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Abstract 

Length scales probed by the large scale structure surveys are becoming closer and closer to the 
horizon scale. Further, it has been recently understood that non-Gaussianity in the initial conditions 
could show up in a scale dependence of the bias of galaxies at the largest possible distances. It is 
therefore important to take General Relativistic effects into account. Here we provide a General 
Relativistic generalization of the bias that is valid both for Gaussian and for non-Gaussian initial 
conditions. The collapse of objects happens on very small scales, while long-wavelength modes are 
always in the quasi linear regime. Around every small collapsing region, it is therefore possible 
to find a reference frame that is valid for arbitrary times and where the space time is almost flat: 
the Fermi frame. Here the Newtonian approximation is applicable and the equations of motion are 
the ones of the standard A^-body codes. The effects of long- wavelength modes are encoded in the 
mapping from the cosmological frame to the local Fermi frame. At the level of the linear bias, the 
effect of the long-wavelength modes on the dynamics of the short scales is all encoded in the local 
curvature of the Universe, which allows us to define a General Relativistic generalization of the bias 
in the standard Newtonian setting. We show that the bias due to this effect goes to zero as the 
square of the ratio between the physical wavenumber and the Hubble scale for modes longer than 
the horizon, confirming the intuitive picture that modes longer than the horizon do not have any 
dynamical effect. On the other hand, the bias due to non-Gaussianities does not need to vanish 
for modes longer than the Hubble scale, and for non-Gaussianities of the local kind it goes to a 
constant. As a further application of our setup, we show that it is not necessary to perform large 
A-body simulations to extract information about long-wavelength modes: iV-body simulations can 
be done on small scales and long-wavelength modes are encoded simply by adding curvature to the 
simulation, as well as rescaling the time and the scale. 



1 Introduction and Summary 



Large Scale Structure (LSS) surveys are becoming larger and larger, and soon they will be able 
to probe cosmological modes whose length scale is comparable to the Hubble scale. General 
Relativistic effects scale as the ratio of the physical wavenumber k/a and the Hubble scale 



and it is therefore important to take these effects into account in order to be able to interpret 
next generation of LSS data. All the relativistic effects are basically projection effects relating 
what happens in one place to what we see: they include such things as lensing, redshift, 
distortion, gravitational redshift, etc. A consistent derivation of them for dark matter has 
been recently performed in pp. Unfortunately we do not observe dark matter directly, but just 
luminous objects. From the observation of them we are able to reconstruct the dark matter 
density field by the realization that collapsed objects are biased tracers of the dark matter 
field. The concept of bias has so far always been defined using the Newtonian approximation 
that is valid for small length scales. The purpose of this paper is to provide a generalization 
of this concept that is valid at arbitrary long-wavelengths. 

Another reason that motivates us to provide such a generalization is due to the recent 
observation that non-Gaussianity in the primordial density field can induce a scale dependence 
in the bias at large wavelengths [21 E]. In the presence of non-Gaussianities of the local kind, 
the bias receives a scale dependence that in the Newtonian treatment behaves as 



where 5 ng is the perturbation to the density of objects, 5 m the perturbation to the matter 
density, and k is the wavenumber of the mode, and where we have neglected factors of order 
unity and the transfer function for simplicity. The important point of this expression is that 
in the presence of non-Gaussianities that have a non- vanishing squeezed limit, such as the 
ones of the local kind or the new ones that have been found in the Effective Field Theory of 
Multifield Inflation jl] with support both on equilateral and squeezed configurations, the bias 
receives a scale dependence at large scales proportional to /nl- This provides an ideal setup for 
measuring non-Gaussianities in LSS, as the signal is peaked on large scales, where theoretical 
predictions are under better control. Indeed current limits on /j^' obtained from the Sloan 
Digital Sky Survey (SDSS) data are already competitive with the ones from WMAP [3], and 
analysis of the bispectrum is expected to be even more promising [5]. 
An odd feature of ^ is that 



It is equally strange that the standard Gaussian bias does not go to zero as k — > 0: one might 
indeed expect that modes longer than the Hubble scale should have no effects on the local 
dynamics. Of course, all of these results are due to the fact that we are trusting ^ way into 
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a regime where it does not apply: as k/a becomes close to H, a proper General Relativistic 
treatment becomes necessary. 

The main purpose of this paper is to provide such a General Relativistic generalization 
of the bias that is valid both in the case of primordial Gaussian and non-Gaussian initial 
conditions. In doing this, we will also provide a way to understand small iV-body simulations 
in the General Relativistic setting, and to show that in order to study the effects of long- 
wavelength modes, it is not necessary to run large, time consuming, iV-body simulations. Let 
us briefly summarize the logic and the main results. 

• Cosmological perturbations become non-linear and lead to collapse only on very small 
scales, where the Newtonian approximation is valid. This suggests that if we insist on 
describing length scales much smaller than the Hubble scale, then the current Newtonian 
description is valid. 

• Given a perturbed Friedman Robertson Walker (FRW) Universe with fluctuations of 
arbitrary length scale, it is possible to identify a coordinate frame valid on spatial 
distances much smaller than the horizon and for an arbitrary amount of time, where the 
metric appears locally as the one of Minkowski space, with small perturbations of order 
(Hx) 2 , x being the spatial distance from the origin. These coordinates represent the 
inertial frame of a free falling observer, and they are called Fermi coordinates [6] . In the 
case where the matter is non-relativistic, in this frame the Newtonian approximation 
is manifest, and we argue that this is the frame where results of small-box iV-body 
simulations can be interpreted. We explicitly construct such a reference frame at linear 
order in the long scale fluctuations for a spherically symmetric configuration of the long- 
wavelength modes, as this is sufficient for the description of linear bias. Generalizations 
to different configurations for the long-wavelength modes or to the non-linear level 
should be straightforward. 

• In these coordinates, all the effect of the long-wavelength mode is included in the map- 
ping from the global frame to the Fermi frame, and in the long-wavelength curvature of 
the local patch. Since for the linear bias we can use spherical symmetry for the long- 
wavelength modes, the long-wavelength part of the Fermi metric must be equivalent 
to that of a curved FRW Universe, and therefore all the effect that a long-wavelength 
mode has on the local dynamics is indeed in the curvature of the local FRW Universe. 
This is given by 

V 2 ((x L ,t L ) 

~ a2R2 , (4) 

where ( is the curvature perturbation in comoving gauge, and here for simplicity we 
have omitted numerical factors given later in the text. 

• This allows us to generalize the concept of bias to the General Relativistic setting, by 
declaring it to be the derivative of the proper number density of objects at a fixed proper 
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time with respect to the curvature of the local Universe: 

, 1 dn p r , V 2 C 



n p dn K np a?H 2 

where n p is the proper number density of objects, 5 np their relative overdensity and 
the dots stand for additional terms coming from various projection effects that we will 
discuss in the text. Here we have neglected numerical factors. This expression makes 
sense physically, as for modes much longer than the Hubble scale, VLk 0, making 
explicit the General Relativistic statement that metric modes that have no measurable 
gradients do not affect the local dynamics. 

In presence of primordial non-Gaussianities, the initial conditions for the fluctuations 
in the Fermi patch can depend on other parameters. In the case of non-Gaussianities 
of the local kind, initial conditions depend explicitly on (, a quantity that has no effect 
on the local dynamics. In this case, we extend the definition of the bias to include the 
derivative of the proper number density of objects with respect to the parameter itself. 
For example, in the case of non-Gaussianities of the local kind, we have 



«, 1 drip V 2 C 



" K-b—^ + b f ^.C + ..., (6) 



where again the . . . stand for additional terms coming from various projection effects 
which we will discuss in the text. We see that the relative factor of k 2 between the 
standard bias and the one induced by /nl is preserved in the General Relativistic limit. 
However, most importantly, the physical effect of long-wavelength fluctuations on the 
local overdensity does not blow up as k — > 0: it is simply the fact that the standard 
Gaussian effect goes to zero while the non-Gaussian one stays constant. 

• Finally we point out that our construction of the local Fermi coordinates shows that it 
is not strictly required to run time-consuming large-box iV-body simulations to study 
the effect of long-wavelength fluctuations: their effect can be simply included by run- 
ning small-box iV-body simulations with different cosmological parameters than in the 
standard cosmology. 

Related works on the way to include long-wavelength perturbations inside small-box N- 
body simulations have appeared in [7J EJ |9]. Related work on the way to derive the bias 
of the local form in the General Relativistic context has appeared in [TOj El H21 [13] and a 
connection of the latter to primordial non-Gaussianities of the local form has been made in 

2 Fermi Coordinates for Perturbed FRW 

Given a sufficiently smooth spacetime, it is possible to identify a set of coordinates cen- 
tered around a timelike geodesic, known as Fermi coordinates |6j. They have two important 
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properties: the metric is approximately that of Minkowski space, with corrections that start 
quadratically in the (space-like) geodesic distance from the time-like geodesic taken as the 
origin, and they are valid in the (spatial) vicinity of the time-like geodesic for all times. 

In an FRW spacetime the Hubble expansion appears in the Fermi coordinates as a small 
correction to the standard dynamics in Minkowski space. This set of coordinates was found 
for unperturbed FRW first in |46j. Here we are going to provide such a set of coordinates 
for a linearly perturbed FRW Universe. We will then argue that in this set of coordinates 
the Newtonian approximation is valid, and that this is actually the frame in which iV-body 
simulations are performed. Furthermore, we will provide a mapping from the local Fermi 
coordinates to the global coordinates of a perturbed FRW, and we will show how simulations 
have to be performed in order to include the effect of perturbations with wavelengths larger 
than the box size. 

Let us therefore find these coordinates. Let us suppose we have an FRW metric with some 
linear long-wavelength fluctuations. We start from a perturbed FRW metric in Newtonian 
gauge: 

ds 2 = - (1 + 2$(£ G , t G )) dt% + a(t G ) 2 (1 - 2#(x G , t G )) dx 2 G . (7) 

In app. [A] we perform the same construction starting from (-gauge. Here the subscript G 
stands for Global to stress that these coordinates are valid for the entire FRW space. A 
great simplification comes from the fact that we wish to study the properties of the large 
scale structures mainly in the regime where the long-wavelength modes are linear: in other 
words, we are mainly interested in the two-point function of large-scale fluctuations. This 
has two consequences. First, the behavior of $ and ^ can be found by solving the linear 
Einstein equations and the linearized equations of motion for matter. For example, we can 
assume that there is no anisotropic stress at linear level, so that ^ = $. Second, if we wish 
to compute scalar quantities (as we will wish), we can use superposition principle to restrict 
ourselves to consider configurations where $ is spherically symmetric around one point, let 
us say the point x G = 0. Generalization to the non-linear treatment of $ is conceptually 
straightforward, but computationally not so, and we leave it to future work [j] 

In order to find the Fermi coordinates (fig. [I]), we can restrict ourselves to the neighborhood 
of a time-like geodesic. Spherical symmetry suggests to consider the geodesic x G {t G ) = 0. If 
we consider modes whose wavelength is much larger than the region of interest, we can Taylor 
expand the metric around the origin, and keep only the leading two derivatives. Notice 
that numerical simulations have to follow dark matter particles, and therefore their region 
of interest corresponds to scales corresponding to the length traveled by the particles, of the 
order of the non-linear scale. We obtain: 

ds 2 ~ - (l + 2$(0,t G ) + $(0,t G ), rGrG r G ) dt 2 G +a{t G f (l - 2$(0,t G ) - $(0, £ G ),r G r G r G ) dx 2 G , 



■""Of course such a non- linear treatment would become much more pressing if we had convincing evidence 
that the primordial perturbation were non-Gaussian. There is some reason of possible excitement: in the 
CMB Gaussianity is excluded only at the 2a level [15] through the analysis of the three-point function of the 
orthogonal kind parametrized by f^ h ° K ' |16j . 
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Figure 1: Fermi Coordinates. 
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We can find the coordinates in which the above metric 



appears in the Fermi way in a simple, but brute force, way that we describe here. A more 
geometric derivation is presented in app. [Bj Let us first warm up by considering the case of 
an unperturbed, curved FRW Universe, whose metric is of the form 



ds 2 



—dig + a(tc) 



cIxq 



[l + \Kxl\ 



(9) 



We consider the curved case here because it will be useful for later purposes. It is easy to 
check that upon the following change of coordinates, valid at small distances 



t G = t L - -H(t L y L 



i 



(10) 



where r\ 



x 



L.I 



+ x \ 2 + x \ 3 an d the subscript l reminds us that these are the Locally valid 



coordinates, the metric takes the form 



ds 2 



l-[H(t L )+H(t L ) 



dt\ + 



H{t L y + 



K 



a(ti 



dx\ 



(11) 



As we had anticipated, for an indefinite amount of time, the metric near the spatial origin 
is approximately the Minkowski one, with corrections starting at order r\ and suppressed by 
powers of Hr^ <^ 1. So for example this metric is valid for distances smaller than Hubble, 
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but it clearly can include cosmologically interesting length scales such as the non-linear scale 
where structures form. 

To consider now the generic perturbed FRW flat space, let us generalize the change of 
coordinates as 

t G = t L - l -H{t L yi - ^ $(o, t')df + 9l (t L y L , (12) 



X) 



X 



G 



a(ti 



l + \H(t L ) 2 rl + f 1 (t L ) + f 2 (t L )rl 



and let us determine the functions /i )2 , <?i, meant to be first order in the metric fluctuations, 
by imposing that the metric in the local coordinates is of the Fermi form, with the additional 
constraint that the spatial part be proportional to 5ij. Notice that we have made the educated 
guess that at the origin the Local time equals the proper time. We will verify shortly that 
this is a good guess. After some straightforward algebra, we obtain 

r tL 

t G = t L - $(0,t L )dt' 



x. 



X 



G 



a(t f ) 



±H{t L ) - H{t L )*$,t L ) - \^t L ), tL - ^ 

t L 

i+$(<5,t L )+H(t L ) I $(d,t')dt'+ 



J%(6,t')dArl, 



{H{t L f + H{t L ) (H{t L f - 2H(t L j) J*J $>(6,t')dt' - H{t L ) 2 ^,t L ) - 2H(t L )®(0,t L )^ r 



.2 
L 



(13) 



Let us recall the common definition of the comoving-gauge curvature perturbation ( 



C(xG,t G ) = -$(x G ,t G ) + 



H(t c 



H(t G ) 



$(x G ,t G ) + 



$(x G ,t G ) 
H(t G ) 



(14) 



and the fact that this is constant for adiabatic fluctuations and for wavelengths longer than 
the sound horizon: 



C(xG,t G ) 



H(t) 



0, (15) 



where the dot stays for derivative with respect to the time variable. This implies that we can 
write ( as 



m = -*(t) - H(t) / dt'$(t') , 



$(t) + H(t)$(t) + H(t) / dt'<Z>(t') = , (16) 



and therefore we can simplify the former expressions to get 
tc = t L — 



X 



G 



a(t f ) 



$(d,t L )dt' - l -H(t L ) (l-$(0,t L ))ri , 

-C(Ml)) • 



;i7) 



The resulting metric is of the form 



ds 2 



1 - H(t L ) + H{t L ) 2 - 2 (H{t L ) 2 + H(t L ) $(0, t L ) - 3H(t L )$(6, t L ) >tL 



-®(0,t L ) itLtL - [2H(t L )H(t L )+H(t L )) I <5>(0,t')dt' ~ 

'0 a \tL) 



dtl + 



+ 



1 - { ~ H(t L ) 2 $(d, t L ) - H(t L )<f>(0, t L ) >tL 



-H(t L )H(t L ) I L $(0,t')^' + 
'o 



r G r G 



a(t L ) 



(18) 



which is valid without assuming that ( is constant. If we use that ( is indeed constant outside 
of the sound horizon, the metric simplifies to 



ds 2 



1 - H{t L ) + H{t L ) 



rc^G 



a(t L ) 



dt 2 L 



(19) 



+ 



H(t L ) 2 , <£>(0,t L ) 

^G^G I J2 



+ 



a(t L y 



dx\ 



The above metric represents the description of a perturbed FRW Universe on scales much 
smaller than the typical length scale of the perturbations. For this reason, in the presence 
of adiabatic perturbations whose wavelength is longer than the sound horizon, it has to be 
equivalent to the local version of an FRW metric, as represented in the local coordinates 



of (j 1 1 p . This is indeed due to Birkhoff theorem. This is in fact true: upon identification of 
e local expansion r 

H L (t L ) = H(t L ) + 



an effective local expansion rate Hl^l) and of an effective curvature Kl given by 

1 



(20) 



Ho,t L ) 



H(t-^ 2 



H(t 



ir '*(ff,t x ) + *<°'**>* 



H{t L ) 



where Hl^l) = ol^l)/^^), the metric (18) takes the form of the curved unperturbed 
FRW Universe in (11) with the simple replacement a — > cll, K — > Kl [^J In this case, the 
local curvature K L is proportional to the Laplacian of the curvature perturbation usually 



'It might be useful to notice that the curvature perturbation K = — £ can be expressed in terms of 



the matter density perturbation in comoving gauge 8 



(com) 
I 



as 



k = [ n ro , + -f ) h$s% 



2 r(com) 



(21) 



where Hq is the Hubble parameter at the present time, £l mj o is the fraction of energy in matter at present 
time, and / = with D being the growth factor such that 5^ com \t) = D(t)5 { com) .The subscript is used 

for quantities evaluated at redshift zero. 
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denoted by £, and is thus constant in time. Hi follows the normal Friedmann equations for 
a curved FRW. 

In summary, we have been able to see that an FRW Universe with a linear adiabatic 
perturbation whose wavelength is longer than the sound horizon can be described, locally, by 
a metric that is very close to the Minkowski one, and is actually equivalent to one of a curved 
FRW Universe. The assumption of adiabaticity and that the wavelength of the mode is longer 
than the sound horizon is necessary in order for the curvature of the Universe to be constant 
in time: it is only in this case that there is one single local history for the Universe, which 
implies that the long-wavelength mode at linear level can be completely re-absorbed into the 
curvature of a local FRW Universe. In practice, this implies that our method of dealing with 
long wavelength perturbations is applicable to adiabatic long-wavelength fluctuations in the 
case where the speed of sound of the fluctuations is very small. This includes a Universe filled 
with dark matter and a cosmological constant, or with dark matter and quintessence with a 
very small speed of sound as the models studied in [17], while it does not apply to models 
with quintessence with non-vanishingly small speed of sound. 

Although for some questions one can restrict to the case of spherical symmetry, this is not 
possible in general. The Fermi coordinates exist also in the absence of spherical symmetry. 
In app. [C]we present the form of the Fermi coordinates starting in Newtonian gauge with a 
plane wave perturbation. 



It is worth to show explicitly how our procedure works in a practical example, where the 
long- wavelength fluctuation is short enough to allow for a Newtonian treatment. Since we 
just said that the effect of a long mode can be re-absorbed in a curvature of the background 
(at linear level and after using superposition principle), this suggests that we should be 
able to re-derive the growth function at second order for short wavelength fluctuations in 
the presence of longer, spherically symmetric fluctuations as derived in [39J in the standard 
perturbation theory approach. Here instead we derive it from the growth of modes in a 
curved Universe. Working only in the limit where all the modes are describable within the 
Newtonian approximation and working in Einstein-de-Sitter space, this calculation is carried 
out in app. [DJ and here we summarize the main results. 

The evolution of short modes in the effective curved Universe is related to the short modes 
in an flat Universe 5 Si fl at as 



where we have restricted ourselves to the matter only Einstein-de-Sitter Universe. In Fourier 
space we get 



2.1 A Simple Check 




(22) 




(23) 



where we have assumed that the long mode is peaked at one frequency ki. 
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Exactly the same expression can be computed in standard perturbation theory as the 
coupling between an general short and a spherically symmetric long mode 5(ki) leading to 

€\ks) + / ^F 2 (k s , k^fyS^ih) = (l + jj^fo)) . (24) 

We see that for long modes sufficiently far within the horizon so that a Newtonian treatment 
is possible, the two expressions agree. 

3 The Coordinate Frame of TV-body Simulations 

Usually, iV-body simulations are performed on very small scales compared to the Hubble 
scale, and no hint is usually given onto in what gauge the calculation is actually performed. 
Further, the equations that are solved in the simulations are not even the General Rela- 
tivistic equations, but the Newton's equations, where all the General Relativistic effects are 
nee lected 

Of course, there is a good reason for this. Usually simulations are performed in boxes 
which are much smaller than the Hubble scale. Since all General Relativistic effects, from the 
corrections to Newton's equation to the specification of the coordinate frame, scale propor- 
tionally to (Ha)/k, these effects are usually negligible. We begin to need to worry when the 
box size of the simulations becomes larger and larger, and reaches the Hubble scale. At this 
point, at least naively, we have to modify our codes to include the General Relativistic equa- 
tions, choose some gauge in which to perform the calculation, take care of what is actually 
the observable quantity that needs to be computed. This is in fact different from Sp m /p m , 
p m being the matter density, as recently stressed in [1] , due to lensing and redshift distortion 
effects. But doing all of this may seem a bit too much: at the end of day, we know that 
large scales evolve linearly, and it is only scales much smaller than the horizon that become 
non-linear and require iV-body simulations. Further, if the sound horizon is much smaller 
than the Hubble scale, local dynamics does not really probe long distances, but it only probes 
distances of the order of the mean free path of the particles, which is the non-linear scale, 
and so it should not be affected by General Relativistic effects. On small scales, we should be 
able to apply the Newtonian approximation, and so our way of doing simulations should be 
fine to describe the small scale non-linearities. There seems to be a tension between including 
long wavelength fluctuations in the simulations, and the fact that the non-linearities occur 
just on small scales. 

This tension has been solved in a recent paper [18] . where it was shown how, exploiting the 
above facts, it is possible to re-interpret the results of current Newtonian iV-body simulations 
directly in the General Relativistic context, by providing a mapping between the results of N- 
body simulations and the fluctuations in a specific gauge valid at arbitrary length scales. The 

3 The fact that usually the spatial coordinates are rescaled by a time-dependent factor equal to the scale 
factor, using the so-called comoving coordinates, should not be misleading: that is just a convenient change 
of variables for the same equations, which are still just the Newtonian ones. 
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only mistakes in this procedure are suppressed by powers of (v/c) <C 1, with no corrections 
of the form (Ha)/k. 

We are now going to argue that this same tension between iV-body simulations and General 
Relativistic effects can be resolved in yet another way, by simply stating that in order to 
include long wavelength modes into the simulations, it is not necessary to make large-box 
simulations, but it is simply necessary to perform small-box simulations, in slightly curved 
backgrounds. The results obtained from the small scale simulations can then be reinterpreted 
as results obtained in local patches of the whole Universe. We will provide such a mapping 

Let us see how this works by showing that simulations can be interpreted in the General 
Relativistic context as nothing but solving the Einstein equations in the frame defined by the 



local coordinates (13), where the metric has the form (11) with the scale factor, the Hubble 



rate and the curvature as given by (20). If we now add short scale perturbations 5$ to the 
metric we have: 



ds 2 



+ 



1 - (H L (t L ) + H L (t L ) 2 j r\ + 25$(x L ,t L 
I / , K 



dt\ 



(25) 



H L (t L y + 



a L (t L ) 2 



25$(x L ,t L ) 



dec r 



As we write down the Einstein equations for a Universe of dark matter particles plus a cos- 
mological constant, where the perturbations are non-relativistic, we immediately realize that 
in the above metric the Newtonian approximation is valid: the metric looks like Minkowski 
with just small corrections, and the system is non-relativistic. Straightforward algebra then 
shows that the Einstein equations take the form of the simple Poisson equation 



V 2 5$(f L ,t L ) = 4nG6p(x L ,t L ) 



(26) 



while the geodesic equation for the dark matter particles takes the form 



6x(t L ) + 2H L (t L )6x(t L ) = -V8$(x(t L ),t L ) 



(27) 



In obtaining the above two equations, we have done several approximations and definitions 
that require explanation. We have defined 



Sp(x L ,t L ) = p(x L ,t L ) 



8nG 



H L {t L f + 



K 



a L (ti 



x(t L ) = H L (t L )x + Sx(t L ) 



(28) 

where here we decided to focus on a ACDM Universe, thought we stress that trivial gener- 
alization of our formulas apply to the case of clustering dark energy [T7j. Notice that the 
unperturbed velocity is nothing but the Hubble flow as seen at small distances from the origin. 



The statement that in order to include large scale modes into small-box simulations one should include 
curvature and a rescaling of the coordinates has been already given in |21j and then more properly in |22j . 
However, a mapping from the frame of the simulations to the global frame had not been given, nor, it seems 
to us, a clear derivation has been presented. Further, all the statements in [2"Tl |2"2"| are not in the General 
Relativistic context. All of this becomes important if we are dealing with modes comparable to Hubble size. 



11 



Then, we have expanded in perturbations by applying the Newtonian approximation: i.e., we 
have counted the perturbations in powers of 5$ ~ v 2 , where v = x(ti), and taken the linear 
equations in these perturbations. Notice that this amounts to taking the leading terms also 
in r\ in the Einstein equations, while we have not expanded in 5p/p. The fact that these 
approximations are justified can be checked a-posteriori, but will become clear in the next 
paragraph. 



In fact, eqs. (26) and (27) are exactly the same equations that are solved in iV-body numer- 
ical simulations. This tells us two important things. First, that the Newtonian approximation 
is indeed justified. Second, most importantly, we now know how to interpret the above equa- 
tions in a General Relativistic setting: they are the equations for a local patch described by 



the local frame. Thanks to the change of coordinates in (13), we can interpret the results of 
the iV-body simulations as points in the full manifold of the spacetime (let us say for example 
as described in standard Newtonian gauge). 

The presence of a long-wavelength mode affects the result of the iV-body simulations in 
two different ways: first it affects the mapping from the global to the local coordinates in 



(13), second it affects the evolution of the short modes by adding a small curvature (20) to 
the effective local FRW Universe. 

In summary, what we found can be synthesized by stating the following simple procedure 
for performing iV-body simulations that include large scale fluctuations. Simulations are to 
be thought of as computing the gravitational structures in the local frame defined by the 



change of coordinates (13). In the presence of a long- wavelength mode, simulations should 
be performed in a curved (background) Universe where the curvature is given by ( |20| ) [^j Any 
scalar quantity measured in the simulations, let us say the proper number density of halos of 
a given mass, should be intepreted as given at this time: 

N— body Simulations — > n,p(xi, Qk{0) > (29) 

where the explicit dependence on ( comes from the curvature, and the superscript L reminds 
us that the output of the A^-body simulations is to be interpreted as given in Local coordinates. 



From the mapping (13), we then finally get the value in the set of coordinates that are globally 



valid, for example in Newtonian gauge: 

np(x G ,t G ;C) = (x L (x G ,t G ),t L (x G ,t G );tt K (()) , (30) 

where the superscript G reminds us that this quantity is defined in global coordinates valid 
everywhere, and we have used that the proper number density is a scalar. 

Finally, we should comment on the initial conditions for the patches corresponding to 
the regions of space simulated in the A^-body simulations. In the case of Gaussian initial 
conditions (we will comment on non-Gaussian initial conditions in the next section), it will 
turn out that to a very good approximation the initial power spectrum, expressed in terms of 
the local coordinates, should be the same as it would be in the absence of the long-wavelength 



5 As we stressed, the same approach can be generalized to include perturbations at non-linear level and to 
compute non-scalar quantities: in this case the local patch will not evolve as a curved FRW. 
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mode. In order to understand the reason of this, it is useful to express the global metric in 
the comoving (£) gauge which is comoving with the density perturbations (see appendix [F]). 
In this gauge, for adiabatic initial conditions, and for modes that are far outside of the sound 
horizon, the metric takes the form 

ds 2 = -dt 2 + a 2 e K dx 2 . (31) 

Let us decompose the fluctuation ( in a long-wavelength and a short-wavelength component 
+ Cs; where i stays for long, and s stays for short. Let us assume for the moment that 
the long-component is on scales longer than the sound-horizon. This means that it entered 
the Hubble scale after matter-radiation equality. In this case, d is constant in time. The 
property of the exponential is such that Exp(C) = Exp(^) + Exp(^), which implies that in 
the limit in which we can neglect completely the gradients of Q, Q can be re-absorbed in a 
constant rescaling of the scale factor, and is therefore unobservable. This implies that the 
local physics (from matter radiation equality to recombination and so on) happens in exactly 
the same way as if the long mode was absent. As we learned in the former section, when we 
consider gradients of the leading effect of the long mode is to induce a curvature for the 
local Universe, which clearly affects the local evolution. So, the initial power spectrum of the 
short scales modes is the one that is obtained in a curved FRW Universe where the curvature 



is given by the Laplacian of Q as in eq. (20). In practice, this means that we should run 



numerical codes as CMBFAST [23J or CAMB [21], run them with the relevant curvature of 
the Universe, and, after a rescaling by the scale factor, simply interpret the output as in local 
coordinates. In reality, it is not even necessary to obtain the power spectrum in such a curved 
Universe as it is easy to realize that the initial curvature is negligible. The relevance of the 
curvature scales as VqQ / '(clqH 2 ) oc 1/af,, and therefore it becomes irrelevant in the past. In 
practice, neglecting the effect of the initial curvature amounts to neglecting terms of order 
Vq0(^,wi)) where tL,in is the initial time of the iV-body simulation. When we later define 
the bias we will define it as the coefficient of proportionality between the local number density 
and V 2 G CiitL, obs), where tL l(} bs is the time of observation. The effect of the initial term scales 
as aL(tL,obs) 2 / ciLiiL,in) 2 and gives a negligible contribution to the bias if the initial time of 
the iV-body simulation is early enough. In practice, this is the simple fact that the curvature 
is irrelevant at early times. This implies that, for long modes that entered the horizon during 
matter domination, the initial conditions for the simulations are equivalent to the ones in an 
unperturbed Universe. 

The situation becomes slightly more complicated for long wavelength modes that enter 
the horizon during radiation domination. In this case, there is a window of time from horizon 
re-entry to matter-radiation equality during which Q depends on time. This means that 
the mode in this case can not simply be interpreted as a rescaling of a and an additional 
curvature term. In this case gradients of the long fluctuation are relevant, as the mode travels 
approximately an Hubble horizon in an Hubble time. In order to evaluate the effect of the 
long mode on the short scale power, one should then solve the non-linear equations that 
couple ( s and Q, along the line of what done in [19]. However, we can argue that this effect 
is negligible. The biasing of structures as due to a long wavelength mode is an intrinsically 
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non-linear effect, and it therefore receives most of its contribution from late times, as density 
perturbations become closer and closer to being non-linear. In perturbation theory, it is 
straightforward to realize that neglecting the non-Gaussianities of the initial conditions set 
up at a time parameterized by a« n amounts to neglecting a non-Guassianity of the matter 
fields at a late time parametrized by a Q b s that is of the order of a in /a i, s . This is equivalent 
to the order of the relative error in the bias we have if we neglect the non-Gaussianity in the 
initial conditions. By taking the initial conditions to be early enough, we can make this error 
small enough. Given the fact that it is quite hard to measure the bias to great precision, the 
initial condition can be set up at a reasonably late time. 

Let us summarize the discussion about the initial conditions. Concerning modes that 
entered the horizon during matter domination, one can simply take the power spectrum in 
local coordinates as in an unperturbed FRW Universe. Concerning modes that entered the 
horizon during radiation domination, one should take non-Gaussian initial conditions that can 
be estimated in perturbation theory as for example in [19]; however, their effect is likely to 
be negligible. The procedure we have outlined in this section enables to extract information 
about very long wavelength modes without practically modifying the iV-body codes, and 
without having to run very large and time-consuming simulations. This should give a valid 
description for certain questions, such as the halo mass function, where spherical symmetry 
that we assumed to derive Fermi coordinates is likely to be valid (see app. [C] for a plane wave 



case). In app. B.3 and B.4 we give a detailed recipe for how to run a iV-body simulation given 
the cosmological parameters and the amplitude of the long-wavelength mode. 



4 Bias in General Relativity and its Scale Dependence 

As an application of our technique we will derive an expression for the bias that is valid in 
the General Relativistic setting. As it has been recently noted in [21 13] in the case of the local 
kind of non-Gaussianities parametrized by the parameter /^l', the bias on large scales (as 
usually measured with respect to to the local matter overdensity) receives a contribution that 
is scale dependent, proportional to 1/k 2 , where k is the wavenumber of the long-wavelength 
mode, proportional to /nl- The same is expected to be true for the new non-Gaussian shapes 
that have been found in the Effective Theory of Multifield Inflation |4J (a generalization of 
the Effective Field Theory of inflation [20J) that have support both in the equilateral and in 
the squeezed limit. These results were derived in the Newtonian approximation, and here we 
will derive their generalization for wavelengths comparable or longer than the horizon. 

4.1 Gaussian Bias 

If we consider surveys that are comparable to the horizon scale, then relativistic effects become 
important and one needs to be very careful in defining observables. We do not directly observe 
the proper number of galaxies at a given point n p (ta,XG) because the photons are deflected 
and redshifted on their way from the source galaxy to the observer. 
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What we can do, is count the number of galaxies in bins of angle and redshift. We will 
refer to the observed number density of galaxies, i.e., the number of galaxies divided by the 
observed volume, as n j, s (z,9,<p). Here z is the observed redshift of the bin, and the tuple 
(9, 4>) represents the observed angular position. 

The observed position (z,9,<p) corresponds to a set of global coordinates (t G ,x G ). Here we 
make use of the fact that a spacetime point can be described in in different coordinate systems 
and that global, local and observed coordinates are just three choices of such a coordinate 
frame that describe the same point. Thus the global coordinates are a function of the observed 
coordinates 

{t G , x G ) = {t G {z, 9, (f>),x G (z, 9, 0)) , (32) 

and since the proper number density n p is a scalar, i.e., a function of the point rather than 
its coordinates, we have 

n p (z,9,(j)) = n p (t G (z,9,(j)),x G (z,9,(j))) . (33) 

To compute the observed number density n obs (z,9,(f>) we need to model both the proper 
density of objects n p and the mapping between proper and observed coordinates. Let us start 
with the proper number density. 



4.1.1 Proper and Observed Number Density 

We have argued that in presence of long wavelength modes, the local inertial frame corre- 
sponds to a homogeneous curved FRW Universe. As a result the proper number density of 
galaxies at the spacetime point is given by the number density in the effective curved Uni- 
verse. We will denote this number n p (t L ;fl K ). The time argument t L stresses the fact that 
the proper time of the free falling observer is in general different from the global coordinate 
time. 

We have: 

n p (z,9,(f)) = n p (t L (t G (z, 9,(f)),x G (z,9,4>))]Q K ) , (34) 

where t L (t G ,x G ) denotes the time in the Fermi frame centered at (t G ,x G ) and VLk is the 
curvature associated with the long wavelength mode. To evaluate this expression we need to 
compute the relation between (z, 9, (p) and ti- We can split this relation in two parts. 

First we can relate (z, 9, <ft) to the global coordinates. As shown in [1] there is a lapse 
between the coordinate redshift 1 + z G = l/a G (t G ) and the observed redshift z 

z - z G = (1 + z G )5z G -> z (35) 

where Sz G ^ z is given in app. [FJ We also need to relate the global time coordinate to the time 



in the Fermi frame at the origin (see eq. 13) 



*l(*g, x G ) =t G + 5t G _> L (t G , x G ). (36) 
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The time shift between the global and the local coordinates is the difference between the 
global coordinate-time and the proper-time. In Newtonian gauge we have 



5t G ^ L (z,e,<f)) =t L (t G (z,e 

1 



t G (z,e, 



$(t G (z,e,<P),x G (z,6,<P))dt' G , (37) 

[C(t G {z,d,<f)),x G {z,e,<p)) +$(t G (z,e,<f>),x G (z,e,<i>))] . 



We can now expand eq. (34) to first order in the perturbations to obtain, 

1 dn Tl 



n p (t L ; £l K ) =n p (t L , VL K = 0) 



1 + 



n p OVLk 



K 



(3f 



where n p is the unperturbed number density at the redshift of observation. Doing so, we have 
performed a split into background and perturbation such that n p (tL, fix — 0) is not a scalar 
but a function of its time argument. Thus 



n p (z, 9, 4>) =n p (t G , tt K = 0) 
=n p (z G ,tt K = 0) 



1 + 



n p dVtx 
1 dn 



Of 



n p dVt K 



p d\ogn p 

"X + TTi 77—, rOZ G ^L 



91og(l + z] 



(39) 
(40) 



where we have rewritten the prefactor and the time shift in terms of the global redshift z G , 
which is possible since there is a one-to-one relationship between redshift and time in the 
auxiliary background Universe that can be translated into a relation between 5t G ^L and 
Sz G ^ L 

x , , a ,s Z G (t L ) - Z G (t G ) 6z G ^ L 

St G ^ L (z,9,4>) = ~ H{z){1 + Z) = ■ ^ 

When calculating spherical averages, the observed redshift z is fixed while coordinate redshift 
z G and global time t G vary. As we will see shortly, it is beneficial to evaluate the prefactor at 

z = z G + Sz G ^ z 



n p (z, 8, 0) = n p (z, Q K = 0) 
We can now define the bias as 



n p o\lK om\\ + z) 
1 dn„ 



(42) 



(43) 



and use that flxit) = 2V G (/(3a 2 H 2 ) in eq. (39). We discuss the relation between this 
definition of the bias and the standard one in the Newtonian approximation in the next 
section. 



4.1.2 Volume Distortion 

Finally, to compute n b s (z,9,(p) we need to take into account the distortions in the volume 
induced by the mapping between (z, 6, <p) and the local frame. These geometric factors were 
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recently derived at linear level in PQ. We denote V p the proper volume corresponding to a bin 
in (z,9,<p) and define 

V P = V P (1 + J), (44) 
where V p is the corresponding volume in an unperturbed Universe and^] 

^ ^ d r l + z r . „ 1 + z dH . 5r , d . 

j =-$-(! + z)—5z G ^ z - 2——5z G ^ z - 5z G _> z -2k + ———5z g _> z + 2— , (46) 
dz tir H dz r 

gives the geometrical projection effects computed in pQ. Finally, we have 



n obs (z, 9, 4>) = n p (z, Q K = 0) 



d lo§ ft 

1 - b nK Q K + p (5z G ^ L - 5z G ->z) + J 

aln(l+2;) 



(47) 



Note that all the terms in the bracket are first order, i.e., they can be evaluated at z, z G or 
t G equivalently, since these agree at zeroth order. 

4.1.3 Observed Overdensity & Averaging 

The observed overdensity is the fractional difference between the overdensity in a certain 
direction and the angular average over the survey area 

d obs {z,9,(f>) = z — . (48) 

n obs [z) 

When evaluating the observed mean number density we can use that all the terms in the 



bracket in eq. (47) vanish, when averaged over a sufficiently big survey area. Hence we obtain 



for the angular average 



f sm9d9d(f> . . . 
n obs [z) = / — n obs {z, 9, (p) = n p (z; il K = 0) . (49) 

J £l SU rvey SUVVey 



Now, the benefit of evaluating prefactor in eq. (47) at the observed redshift becomes obvious 



Since the observed redshift is fixed, n p (z, VLk = 0) agrees with the survey average and we have 
for the observed overdensity (we will ignore the additional effects on monopole and dipole, 
which are influenced by the contributions at the observer's position), 

O lo£[ Tl 

6 obs (z, 9, <P) = -ba K n K + g ln(1+ P z) (Szg^l - 5z G ^ z ) + J . (50) 

The volume distortion is in principle observable and thus has to be gauge invariant by itself. 
The first term — ^q k is the number of collapsed objects in the inertial frame and thus 



6 Our expression for J assumes that the survey is volume limited. If instead the survey is flux limited, we 
have to add the corrections due to the change in the apparent luminosity. In this case we have to replace J 
with 

J -> J-5p6V L . (45) 



See app. |F]for details. 
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totally independent of the choice of coordinates on the global manifold. The remaining 
redshift lapse is gauge invariant as we show in app. [F} Together with the first term it forms 
another observable. With the above results the observed number density can be written as 

n obs (z, 9, 4>) = n obs (z) [1 + 5 obs (z, 9, 0)] . (51) 



We can also relate the expression in (47) to the overdensity in the global coordinates, 5a(ta, %g) 
(n p (tG,xa) — n p (ta)) /n P (tc), where the averaging is done over hypersurfaces of constant co- 
ordinate time. We obtain 



n obs {z,9,4>) = n obs (z) 



' d\ogn p ' 

1 4" Oq ^"j — 7T— — rdz G ^ z + J 
alog(l + z) 



(52) 



Note that in the case where the tracer has a number density that scales like (1 + z) 3 the 
combination 5q — (<91ogra p /<91og(l + z)) Szq^z becomes 5q — 35zg^ z m agreement with pQ. 

In our formalism it was natural to define the bias directly in terms of the Laplacian of 
the ( perturbation at the point of interest, which in turn is proportional to the curvature of 
the local FRW Universe. Because of the Friedmann equation, the curvature turns out to be 
proportional to the overdensity of the Universe at the source galaxy position, as shown next. 
This offers us a procedure to extract the bias from iV-body simulations: run simulations with 
varying Q^, and then take the derivative with respect to this parameter. 



4.2 Comparison with Standard Newtonian Treatment of Bias 

Our bias definition tells us that we should take the derivative of the number density with 
respect to the curvature of the local Universe. While our receipe is well defined in the full 
General Relativistic setup, it still should agree in the limit in which the long mode is well inside 
the horizon, so that the Newtonian approximation is valid for the long mode itself. However, 
in this case a naive look at the expression might make us think that the two procedures do not 
agree. Indeed, in the classical Newtonian treatment, the bias is defined as the derivative of 
the number density with respect to the local long-wavelength overdensity. In this section we 
will first relate the above bias definition to an overdensity and then consider the subhorizon 
limit. 

The curvature energy density of the local Universe scales as Qk = ^k,oHq / (clgHg) 2 and is 
thus fully specified by its value at redshift 0. The latter can be related to the matter density 
in (synchronous) comoving gauge as 

O - K 2V 2 ( _ f f H$ \ , {com) 

" Acorn) 
7,0 

where the last two equalities are valid for a Universe with time varying dark energy and a 



n m ,o + §/oUK °r (54) 



ACDM Universe, respectively. Hence the bias term in eq. (47) can be written as 
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We restored the time dependence of the long wavelength density perturbation, dividing by 
the linear growth factor D(t). Our new bias bn K is related to the standard bias parameter by 
a time dependent but scale independent factor. From the above equation we can see that the 
density perturbation in the comoving gauge is equally suited, at an algebraical level, as an 
expansion parameter for the galaxy bias, but the justification of this statement relies simply 
on the proportionality of b to bu K . Further, the bias expressed in terms of Qk makes manifest 
its gauge-invariant physical origin and the fact that the biasing vanishes for modes longer 
than the Hubble scale. 

Well inside the horizon (k 3> aH) the velocity term in the relation between comoving and 



Newtonian gauge matter overdensity (see eq. (156) in app. |Fj) becomes negligible and thus 
both density perturbations reduce to the Newtonian density perturbation 8f ~ <5 / com ~ 81. 
Furthermore, inside the horizon the volume distortion as well as the lapse between the global, 



local and observed redshift are negligible. Thus eq. (50) reduces to 



6 0ba {z,8,<l>)=-b slK {t)n K (t) = b(t)6 l (t) , (56) 

which is the standard relation between observed tracer overdensity and underlying matter 
overdensity in the Newtonian approximation. 

Finally, we point out that another way to understand the connection between our bias 
bn K and the standard one is by referring to the peak background split method. There, in 
the Newtonian context, it is usually assumed that the presence of a long scale mode can be 
interpreted as a shift of 8 C : 8 C — > S c — 8i, and after Taylor expansion we obtain the expression 
for the linear bias. In our context, the presence of a long mode is instead interpreted as 
a curvature of the background Universe, and therefore we have to rescale 5 C accordingly to 
8 C (VL K = 0) — > 8 c (Vt K 7^ 0) and then Taylor expand. In app. [Ej we show that indeed the two 
approaches are equivalent on short scales. 



4.3 Bias in Presence of non-Gaussianities of the Local Kind 

So far we have assumed that the only way a long-wavelength mode affects the local structure 
formation is through its dynamical effects: that is by changing the local geometry and by 
introducing curvature in the resulting local FRW Universe. If the initial conditions are Gaus- 
sian, this accounts for all the effects of the long mode on local processes: in the linear regime 
the statistical properties of the short wavelength modes are decoupled from long wavelength 
modes, and the non-linearities kick in only at late times on small scales, where all the effect 
of the long mode can be absorbed by a redefinition of the local expansion history. If the 
initial conditions are non-Gaussian, then the statistical properties of the initial short scale 
fluctuations are in general affected by the presence of a long mode and this has to be taken 
into account. The scales that become non-linear are very small compared to the horizon, 
and the scales that we are interested in are much larger than the non-linear scale. Thus, in 
order for the properties of the short scale fluctuations to be affected by the long mode, the 
non-Gaussian initial conditions need to be such that they correlate very long and very short 
modes. 
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In general the description of the statistical distribution of modes in the initial conditions 
requires knowledge of all the moments of their distribution. For special cases a limited set of 
parameters p is sufficient. For instance, if the initial conditions are Gaussian, they are fully 
quantified by their variance. If the parameters p depend on the long wavelength amplitude, 
then the proper number density of objects has an additional explicit dependence on the long 



wavelength amplitude. Thus we can generalize eq. (30) to: 



n p {x G , t G ; C) = n p (x L (x G , t G ),t L (x G , t G ); Q K , p(C)) . (57) 

These parameters p represent all the relevant information needed to describe the initial con- 
ditions on small scales. The abundance of objects of a given mass M is mainly sensitive to 
the amplitude of fluctuations smoothed on a scale enclosing the mass, given in terms of the 
variance ou- There is also a weak dependence on the slope of the power spectrum at the 
scale M and possibly on parameters describing deviations from a Gaussian distribution of the 
small scale modes, e.g. skewness. For definiteness we will consider only the dependence on 

The so-called local kind of non-Gaussianities [25] that can be produced in multifield in- 
flationary models j26j E7] or in the new bouncing cosmology [28J provides an example where 
(7m depends explicitly on the long wavelength amplitude ( Q In these models the initial 
conditions are such that the curvature perturbation is a non-linear function (local-in-space) 
of an auxiliary Gaussian random variable ( g : 

3 



C(*g) = Uxg) ~ {( 9 (x G ) 2 - (Q) 



(5f 



If we decompose ( into long and a short modes as we did before, we can see that the short 
mode takes the form 

C. ^ (i - ^nl Cm) C 9 ,s , (59) 

where we have neglected a term in £| which is irrelevant for our discussion of the bias. From 
this equation we see that the variance of the short scale power is modulated by the long mode. 
This implies that in the case of local non-Gaussianities there is an additional source of bias. If 
we set up the initial conditions for the simulation in the presence of non-Gaussianities of the 
local kind, the resulting proper number density of halos n p will depend on the long-mode not 
only through its explicit dependence on the curvature of the local Universe, but also through 
the dependence on the initial power spectrum of the modes [^J Eq. (47) is generalized to 



n obs (z,9,4>) 



n p (z; n K = 0, p) x 
1 - bo,. 



2V 2 G C , 1 dn p da 2 M 
nK 3a 2 H 2 



n p da 2 M d( 



3 log TL 

C + ol ft , P \ ( Sz G->L ~ Sz G ^ z ) + J 

alog(l + z) 



n p {z) 



d log n p 
<91og(l + z) 



(Sz G 



Sz, 



G- 



+ J 



(60) 



7 The same is expected to be true for the new non-Gaussian shapes that have been found in the Effective 
Theory of Multifield Inflation that have support both in the equilateral and in the squeezed limit [3] . 

8 Though we are now talking about non-Gaussian effects, notice that we are consistently treating the long 
mode at linear level. 
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Figure 2: Observed galaxy power spectrum for z — 1, bn K = 1-5 {b = 2) and <91ogn p /<91og(l + 
2) = 3. We choose the following cosmological parameters: Q m = 0.28 , = 0.84 , H = 0.70. 
Le/f panel: We show the spectra parallel to the line of sight (red) and transverse to the 
line of sight (blue). The solid line is for Gaussian initial conditions, whereas dot-dashed is 
/nl = +0-5 and dashed is fj$£ = —0.5. The lower black line is just the power spectrum 
of density in comoving gauge, the upper is multiplied by the redshift space distortion factor 
(1 + f /b) 2 to give the power parallel to the line of sight. We see that the effects of non- 
Gaussianity and GR-effects on the power spectrum differ, because the latter depend also on 
the line of sight parameter /i through the peculiar velocity effects. Right panel: Same as left, 
but orange lines show non -Gaussian power spectrum without the GR-effects (just redshift 
space distortions). 



where da 2 M /d(, = —12/^ /5 is independent of M and p describes the initial conditions in 
absence of long perturbations. We see that in presence of non-Gaussianities of the local kind 
the bias receives an additional contribution proportional to (, while the standard Gaussian 
contribution is proportional to V 2 C- There is a relative scale dependence proportional to k 2 
between the two. But this does not imply the very unphysical result that the bias blows up as 
k —7- 0. It is rather the fact that the bias for large scales should be interpreted as a different 
bias: as the coefficient of proportionality between the local number density and ( and V 2 C 
The final expression for the observed overdensity in presence of local non-Gaussianities is thus 

5 obs (z,0,<P) = -bn K Q K + kC + g ln(1+ P ^ (Szg^l - Sz G ^ z ) + J . (61) 

The presence of $ terms in the redshift lapse terms and the volume distortion term mimicks 
/nl °f order unity. But this should not bias any measurement of non-Gaussianity since the 
General Relativistic effects are calculable and can thus be removed from the measurement. 

9 Our conclusions about the bias as due to local non-Gaussianities are in general agreement with the ones 
of [12] , though they differ in the way the results are derived and in parts of their interpretation. We stress 
that our derivation does not crucially rely on the assumption of spherical symmetry. It should allow for a 
straightforward generalization to the non-linear case where spherical symmetry can not be used. 
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As in the case of Gaussian initial conditions our formula can be applied directly to the 
results of iV-body simulations, but for illustrative purposes we can also calculate the effect 
analytically assuming that the number density of collapsed objects is described by a universal 
mass function 

n p ^f(^-) , (62) 

i.e., it is a function of the peak height, the ratio of collapse threshold and fluctuation ampli- 
tude. In this case the derivatives of n p with respect to a 2 M and with respect to the curvature 
V 2 C are related: 

drip _ dn p da 2 M l dn p S c da 2 M 

d( ~da 2 M d( ~ 2d5 c a 2 M d( ' 1 j 

dn p drip dS c 



d(-n K ) d6 c d(-n K ) ' 

which means that the Gaussian and the non-Gaussian bias are analytically related. 

For general initial distributions, the additional contribution to the fluctuations in the 
proper number density arise from dp/d(. In the case of <jm considered above, this is nothing 
but the squeezed limit of the three point function, because this measures the coupling between 
small and large scale modes. Squeezed in this context refers to the fact that we are talking 
about a correlation between short and long wavelengths so two of the momenta in the relevant 
three point function are very large compared to the other one and thus the three momenta 
form a squeezed triangle. In the local model that we use as an illustration, the derivative 

da 2 M , , 

(64) 



d((k) 

is independent of both k and of M (we have explicitly pointed out that the derivative might 
be different as a function of the wavenumber of the long momenta). Relatively simple models 
can and have been constructed where this derivative depends on both M and/or k [29]. Even 
when this derivative is not constant our formulas remain valid. 

As an illustration, in fig. [2] we show an example for the the observed galaxy power spec- 
trum for bQ K rj 1.5 (b = 2) at z — 1 assuming a volume limited survey. The plots show 
the power parallel and orthogonal to the line of sight for a sample with evolution slope of 
<91ogfip/<91og(l + z) =3. We are adding local non-Gaussianity of = ±1. The right 
panel shows that ignoring the GR-effects could lead to a fake detection of /nl = Oil), but 
this degeneracy is broken if modes both transverse and along the line of sight are considered. 
This is because GR-effects have a peculiar velocity contribution that has a /i dependence, 
where /i = cos# and 9 is the angle between the Fourier mode angle and the line of sight. 
For a related study on distinguishing GR effects from primordial non-Gaussianity see [30] . 
Note that the magnitude of the GR-effects depends on the redshift distribution of the sample. 



The non-Gaussian bias parameter is calculated from the Gaussian bias using eq. (63). For 
the evaluation we are neglecting all the line of sight integrals (convergence, Shapiro-delay, 
integrated Sachs- Wolfe effect), which contribute power mainly to transverse modes. For the 
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details of the evaluation of the observed power spectrum, we refer the reader to eq. (174) in 
app.[F} 

It is also important to note that the relevant quantity is the change in the amplitude 
of fluctuations at a given physical scale M not of course a comoving scale. In single field 
inflationary models this derivative goes to zero in the squeezed limit, when k corresponds to 
a much larger scale than M. In fact it goes to zero as the square of k just because the long 
wavelength mode affects the production of the short modes during inflation only through tidal 
type effects. In a sense it goes to zero in this way for reasons identical to the ones that lead 
to the VgC dependence in the bias formulas. Thus, in single field inflationary models there is 
no modulation of the proper number density that scales with lower powers of k than V^CQ 
The reader familiar with the standard calculation of the single-field inflationary three point 
function might recall that in the squeezed limit they do not seem to vanish but that they 
satisfy a consistency condition where the shape of the three point function looks like that of 
a local model with an amplitude given by the tilt of the fluctuations usually called (n s — 1). 
But this dependence arises entirely from the fact that what is being calculated is a three 
point function in terms of comoving momenta. If expressed in terms of physical momenta, 
the (n s — 1) is exactly the amplitude required to make the relevant derivative vanish. 

4.4 Observing Local-type non-Gaussianities in the Presence of GR 
Corrections 

The salient fact about the local-type non-Gaussianities is that they induce a dependence of 
the proper number density of objects on the long wavelength modes that is much stronger 
than what the dynamical effects can produce, proportional to ( rather than V^C- Unfor- 
tunately when we count objects in our Universe there are projection type effects that make 
the observed densities depend directly on ( even if the proper density does not. The vol- 
ume corresponding to a given observed range of angles and redshifts varies as a result of 



the long wavelength modes and results in the factor of J in eq. (47). Furthermore a given 
observed redshift corresponds to a different proper time in different directions resulting in 
the terms proportional to <91ogn p /<91og(l + z) . Both of these terms lead to contributions 
proportional to (, contributions that have the same form as that coming from the local-type 
of non-Gaussianities. Failing to correct for them would bias the results for /j^' by a number 
of order one which depends on the details of the population of objects surveyed. 

Of course the various terms have different dependences on the properties of the objects as 
they depend on different derivatives of rip . The effects will also depend differently on redshift 
and furthermore, because the GR effects are projection effects induced by the intervening 
matter, it may be possible to distinguish them using observations of the distribution of matter 
at the intervening redshifts. It is beyond the scope of this paper to quantify the extent to which 
these different effects may be isolated in practice or what it is required of the observations to 
distinguish them. 



10 One can construct examples were there is an intermediate "squeezed regime" over which the scaling is 
different than k 2 but for sufficiently large ratio the scaling needs to be k 2 [3TJ [32, 33, 34| l35| 136] . 
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It is clear however that the GR effects are just projection effects, so if we were able to 
construct observables that were directly sensitive to quantities in the local frame we could 
side track those difficulties. In this section we just want to point out that this is in principle 
possible. We will not address wether this can be done in practice given our current tools 
or wether this route is better than just trying to correct for the projection distortions in a 
realistic situation. 

To be able to ignore the projection effects we would need to be able to measure the proper 
density of some object at a given proper time. Thus we would need a ruler that would allow 
us to measure distances independently of the observed angles and redshifts and we would 
need a clock that would allow us to compare regions of the Universe at the same proper time 
independently of the observed redshift. If we managed to find such local clocks and rulers the 
observed density should only depend on VgC in the absence of primordial non-Gaussianity. 
In fact there should only be a Vq( dependence in any single field model of inflation. 

There are many such rulers that one could imagine using. One option is to use the acoustic 
scale. This could be used for example by measuring the the number of objects in regions of a 
given size in units of the acoustic scale. The acoustic scale can be determined by measuring 
the correlation function of these or other objects. Another option is to measure the ratio of 
the densities of two tracers. Then the volume projection effects would cancel, in a sense we 
are using the density of one of the objects to define the ruler for the other. 

We still need a clock to make sure that one is comparing the number densities at a 
fixed proper time rather than observed redshift. This difference is responsible for the terms 
proportional to (<91ogn p /<91og(l + z)) in eq. (47). This appears a bit more tricky but not a 
problem of principle. One needs to date the object observed independently of their redshift, 
something that happens automatically for tracers that appear only at a characteristic time 
in the history of the Universe. Examples of such things might one day be the first stars or 
perhaps quasars could be used as their abundance has a peak in redshift. In other words, 
the ratio of densities of tracers that come from a given proper time and could be identified 
without using the observed redshift would only depend of the long wavelength modes through 
the V|C- 

A similar construction for measuring the three point function in the squeezed limit could 
be accomplished using the CMB. The CMB comes already from a defined proper time, the 
recombination of hydrogen provides the clock. So one could use the dependence of the small 
scale power on large modes as a test of the squeezed three point function. One should use 
a local definition for fluctuations and power, meaning normalizing the fluctuations to the 
mean fluctuation level in the region of interest to eliminate the equivalent of the Szq^z term 
in our equations for the densities of haloes. There is still the projection effect related to 
the mapping between angles and physical distances at recombination. This however can be 
avoided by comparing the amplitude of fluctuations at a fixed scale measured in units of 
the acoustic scale, thus at a fixed physical scale. For example the amplitude of the power 
spectrum at the iV-th peak should only depend on the long modes through the V^C i n the 
absence of primordial non-Gaussianities. One could also use the anisotropies in the small 
scale power to de-lense the CMB along the lines considered in [371 EE] . 
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Appendix 

A Fermi Coordinates from (-gauge 

Here we give the change of coordinates necessary to go from £-gauge to the Fermi coordinates 
in the case of a spherically symmetric perturbation. 
In £-gauge the metric takes the form 

ds 2 = -N 2 dt 2 + 5ije K a 2 (dx i + N { dt) (dx j + N j dt) , (65) 

where we have used the ADM parametrization. In this gauge time diffeomorphisms are fixed 
by requiring Tf = 0. The lapse A" and shift N l are constrained variables, whose solutions in 
terms of ( are [13] 

at -i I C AT V G ,iC H a 2 V G ,i i (aa , 

n = 1 + h> N * = -^r~ip^vi c - (66) 

The equation of motion for ( reads [13] 

1 n fa 3 H \ H V 2 r ( , , 

Outside the sound horizon and assuming c s constant, we can simplify it to 

t _ H c 2 V 2 G ( 



(^(A)+3f) a2 



Plugging back in (66), we can simplify the expression for A^ and Ni to be 

/ 

K H ' H " SH -d t [k/H 2 



N ~ 1 , Ni 



1 H 1 

+ 777 t. r V ^C • (69) 
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At this point we proceed as in the main text. We Taylor expand £ around the origin assuming 
spherical symmetry, we make an ansatz for the change of coordinates, and we impose the 
resulting metric to be in the Fermi form. After some straightforward algebra, we obtain for 
the change of coordinates 



' G 



X G ~ 

and for the metric 
ds 2 = 



tr - 



X, 



a(t L ) 



H{t L ) + 



C,rGi 



r G 



1 - 



HH 



-3H 4 - 2H 2 + HH 



L 5 



1 + 



H(t L )\ 2 



1 " C(0) 



(70) 



H(t L ) + H(t L 



■ 2 C,r G rG 



1 



a 



H 2 [ZH 4 + 2H 2 - HH 

3 TJ I Tj2 TT I o Ej3 i TT2 \ i Tji 



9H*H + 9H 6 H 2 + AH - 3H 



-QH 5 HH + 2H 3 H 2 H - AHH 3 H + H 2 H -2H 3 + H 2 ) + H 4 12H 3 + H 2 - HH 



+ 1 



H{t L ) 2 C 



H 2 H 



2 a{t L ) 2 -3H A — 2H 2 + HH 1 ' ' "'' ' 



As expected, this metric has the same form as the Fermi patch of a closed FRW Universe 
with 



H L (t L ) = H(t L ) + ^(i-, 



HH 



a(t L ) 2 \H 3H 4 + 2H 2 - HH 



(72) 



K L = --v 2 G ao,t G ) 



in agreement with what found in the Newtonian- gauge case. 



B A Geometric Derivation of the Fermi Coordinates 

In this section we will describe how the Fermi coordinates can be constructed from a geometric 
point of view. 

The starting point for the derivation will be a free falling observer moving along a timelike 
geodesic ft, (7) in the background Universe (fig. [3]). His coordinate axes are described by an 
orthonormal set (e*o, e\, e*2, e^) which is parallely transported along h. Thus if e*o is tangent 
to the geodesic ^(7) at its origin it will remain so for all values of the affine parameter 7. 
Without loss of generality we can assume e to be timelike and the e,, i = 1,2,3 to be 
spacelike, and the geodesic to be the origin of the global coordinate frame x G = 0. 

Now we consider a point P = ft (70) on this geodesic. Our goal is to describe the spacetime 
in a neighborhood U of P starting from the global metric at P. Any point Q in the vicinity 
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h(Y) 



Figure 3: Geometrical Construction of the Fermi Coordinates. 

of P can be connected to P with a geodesic g(X) that is perpendicular to the tangent vector 
of h at P, i.e., its tangent vector v at P is a linear combination of the e*j. The coefficients 
of this linear combination are the Fermi coordinates and the time component of the Fermi 
coordinates is chosen to be the proper time of the observer moving along h. The point Q can 
thus be fully described by the proper time r of the observer at P, the direction cosines x l and 
the length of the geodesic A joining Q with P. For simplicity we will normalize the direction 
cosines such that the point Q corresponds to A = 1. This prescription is the natural extension 
of the flat space polar coordinates to curved space. The observer points in a certain direction 
defined by the direction cosines x l L and then follows the geodesic defined by the direction. 
The initial conditions for the geodesic connecting P and Q can thus be summarized as 

x\\ = 0) = , r(A = 0) = t L , (73) 



die 

dX 



0, 



dx 



G 



A M d ^ 



x l L ei(t Q ) 



A=0 



The point Q G U with Fermi coordinates x^ L is then found by propagating along g(X) un- 
til A = 1. We now have to find the mapping between arbitrary coordinates x M and the Fermi 
coordinates defining the geodesic g(X). This can be done by solving the geodesic equation 
for g(X) 

perturbatively using the power law ansatz 

x"(A) =(4 + o%\ + a£A 2 + c^A 3 + . . . . (75) 

The validity of this series is clearly limited as is the validity of the Fermi coordinates them- 
selves, which is obviously related to the curvature of the spacetime. The four vector formu- 
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lation for the initial conditions stated above is 



off =(*„, 0,0,0), 
dx^ 



(76) 



d\ 



A=0 



where to is the coordinate time corresponding to 70. The coefficients of the second and third 
order terms in the Taylor series follow straightforwardly from the geodesic equation evaluated 
at P 

„ 1 d 2 x» 



a, = 



2! d\ 2 
1 d 3 x» 



A=0 



(77) 



6! d\ 3 



A=0 



dT 11 



where we already simplified using the initial conditions. In the following two subsections we 
will describe the mapping for two specific cases: perturbed and unperturbed FRW Universes. 



B.l FRW 



We will now follow the above procedure for the homogeneous Friedmann-Robertson- Walker 
metric 



ds 2 = -dt 2 G + a(t G ) 



dx 



G 



[1 + \Kx 2 G ] ' 

The vierbein associated to a comoving geodesic is 

[e o r = (1,0,0,0) , 
[e 2 ]" = a- 1 (0, 0,1,0) , [e 3 

The linear coefficients in the geodesic expansion read 

a£ = a(t ) _1 (0,a: L ,|/ L ,2; L ) . 



= a-\0,l, 0,0) , 
" = a _1 (0, 0,0,1) . 



(7? 



(79) 



(80) 



Hence, the first order spatial separation is x G ~ x % L /a(t ) and thus x % L is nothing but the 
physical separation of Q from P. Up to third order in the affine parameter we obtain 



to —tL — 



■X / 



x 



G 



a(t L ) 



Ex\ 



1 + 



(81) 



E 2 x\ 



This leads to the following metric in Fermi Normal coordinates 



ds 2 = - 



+ 



l-{H{t L ) + H 2 {t L )^ 



dt 



Sij - H 2 (t L ) + 



dx j^dob £ 



(82) 
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The above metric has non-zero off-diagonal contributions. The general transformation to 
remove off diagonal terms can be derived considering the metric in the old coordinates x 



and new coordinates x(x) 



ds 2 = A Sij dx l dx^ + B XiXj dx l dx? , 



ds 2 = A dx % dx* 



Using the ansatz x % = x l (l + 7a; 2 ) we obtain the condition valid at second order in x: 



B 



7 = 



AA 

A — A(l + 2 7 x 2 ) 



(83) 
(84) 

(85) 



For the FRW case we have 7 = —H 2 /12 — K/(12a 2 ) and the time component up to second 
order is unaffected 



tc —ti 



H{t L ) 



x 



a{ti 



2 

1 + 



-x 



L 5 



L ,1 , H{tL)\ 2 



-X, 



L > 



finally leading to the following metric 



ds 2 = - 



1 - (H{t L ) + H{t L ) 2 ) 



x 2 



dt 2 L + 



1 - H{t L ) 2 + 



K 



a(t L yj 2 



which has the desired form. 



X I 



dx L , 



(86) 



(87) 



B.2 Perturbed FRW 

Let us now consider a perturbed FRW Universe in Newtonian gauge 

ds 2 = -(l + 2$(t))cft 2 + a 2 (t) (l - 2^{t)^dx 2 . 

We assume vanishing anisotropic stress leading to $ = ^ . The vierbein associated to the 
coordinate frame is 



[e p = (l-$, 0,0,0) , 
[e 2 ]" = a- 1 (0,0,l + $,0) , 



[e 1 ]^ = a- 1 (0,l + $,0,0) , 
[e 3 ] M = ^(0, 0,0,1 + $) . 



(89) 



We can now for simplicity expand the potentials around P 



$(x,t) = $(0,t) + 



19 2 $ 



2dr 2 G 



r 2 G = $(d,t) + h(6,t), rGrG r 2 G , 



(90) 
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where we asumed spherical symmetry [^j leading to 

ds 2 = - (l + 2$(0, t) + $(0, t) jrGrG r 2 G yt 2 G + a 2 {t) (l - 2$(0, t) - $(0, t),r r r&)d^, • (91) 

There is no linear term in this expansion, because we require the potential to be different iable 
at r = 0. Let us proceed to find the Fermi coordinates. As noted above, the Fermi time is 
the proper time of the observer following the central geodesic h 



t L = [ V^dt' = t+ [ <S>(6,t')dt' 
Jo Jo 



(92) 



The coordinate time at P thus is to = t^ — f Q L <J>(0, £') dt', where in the integral boundary 
ti = to at leading order in $. This leads to the following expansion factors 

<=(t , 0,0,0) , (93) 

„ 1 + ^(0, to) fn s 
a(t ) 

At this point, simple algebra as shown in the former section leads to the same relationship 



among the coordinates as in (17) and to the same Fermi metric as in (19). 



B.3 Local Expansion Factor 

With the aim of giving very specific recipe for running simulations given a certain long wave- 
length fluctuation, we provide some more specific relations. Some expressions can be simpli- 
fied by noticing that the potential in Newtonian gauge and the density perturbation in the 
comoving gauge are related by (see Appendix [F]) 

V 2 $(f,t) = ^ixGa 2 p5\ com \x,t) . (94) 

Note that this equation is exact, even on horizon scales. We define the growth factor in 
comoving synchronous gauge as 5\ com \t) = D(t)5l c Q om \ which is normalised to unity at present 



time. We also define the logarithmic growth factor f(a) = d\n D/d\na. From eq. (94), we 



define the growth factor of the Newtonian potential D as follows $(£) = D(t)&o/ a (t), where 



$o is the present day value. From the linear growth and eq. (94) it follows 



®(0,t),tr G ra = #$(0, t),r G r G (/ ~ 1) • (95) 

Using the constancy of ( we have shown that the value of the metric perturbation at the 
origin is irrelevant for the local expansion. Thus it only remains to derive the rescaling of the 



11 Note that 

2 2 9$ <9 2 $ <9 2 $ 

V $ = 1 = 3 

r dr dr 2 dr 2 

where the last equality is true for a power law in r assuming no linear dependence in r. 
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expansion factor corresponding to the effective local Hubble rate. Starting from (20), which, 
by defining H L (t) = H G (t) + 5H(t), gives 

1 H G (t) ( H ®(d,t), t \ 

\ / ,r G r G 

we can find the corresponding rescaling for the expansion factor using the ansatz 

a L (t) = a G (t) (1 + 8a(t) rel ) , (97) 
where 5a re i has to satisfy the following: 

fa "-' f > = 3aM f )g G W V °( $(0 " i) + C(0 " f> ) ' <98) 



Sa re i(t) 



where we have chosen the constant so that the two scale factors agree at early times. 5a(t) re i can 
be numerically integrated from the transfer functions for any given cosmology. 
Finally, the Friedmann equations in the Fermi frame read as 

8ttG_ 8ttG_ K 

A 

and 



H 2 L = ^PL + ^Pde - ^ , (99) 

O O CL j 



a L AtxG 8nG 

— = ^PL + 7, PDE , 100 

where pi is the local mean matter density. From the rescaling between the global and local 
Hubble rate we can derive the rescaling of the local mean density 

HL{tf + ~ *fp DE (t) _ p L (t) 

H G {tf - 8 -fp DE (t) ~ pa(t) ' 1 j 

leading to 

PL® = Pc{t) + ^I'gf = Pa® (1 + *\ COm) (*)) • (102) 

This relation can be intuitively understood in the Newtonian context: the long wavelength 
density just rescales the local mean density. This relationship gets upgraded to the relativistic 
setup by using the comoving gauge overdensity. 

For definiteness we give also the closed form expressions for the expansion and Hubble 
rate in a ACDM background 



i 1 - tmm = Hcit) v - i ms ' (t) 1 • (io3) 



a L (t) =a G (t) 
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B.4 Local Density Parameters 



The time evolution of the local patch is determined by the local Friedmann eqns. (99) 
and (100), which are parametrized by the effective local density parameters. We will now 
provide the explicit mapping from the global to the local cosmological parameters that are 
needed for simulations. We specialize to ACDM for simplicity though, as we stressed, our 
approach applies also to clustering dark energy. From the relationship between Hq and Hl 
and from the definition of of K given in (20), we have 

K 2 1 



n K>L (t L ) 

n A , L (t L ) 



a L (t L ) 2 H L (t L y 3a L (t L ) 2 H L (t L ) 2 



v 2 G C(£ G ) 



(104) 
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n m , L (t L ) = 1 - Q K (t L ) - Q A {t L ) 



It is convenient to normalize the cosmological parameters at a = 1 (for us = 1), which 
leads to 
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(105) 



° n L,0 



m,L,0 



1-Q 



^A.O 



where the subscript o stays for evaluating the quantity when at, — 1. In order to be able to 
use the above formulas, we simply need to find the time at which a& = 1. This can be 
found by solving eq. (97) with a L = 1 to identify t£ )0 . From there, by plugging into (20) we 
get Hifi- The former expressions can be further simplified to give 
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The local Hubble rate at = 1 is given by 

I V 2 G C(0) 



Hl,o = H G)0 1 + 
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(107) 



As an example we consider the WMAP5 Flat ACDM cosmology with matter density pa- 
rameter VL my Q = 0.28 and Hubble constant Hg,o = 70kms~ 1 Mpc" 1 . For a long wavelength 
amplitude of (Jjq 1 = 0-1 corresponding to Vq(/Hq = 0.89 we obtain 



m,L,0 



0.30 
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K,L,0 



-0.06 



n 



A,L,0 
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(108) 
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Figure 4: Time dependence of the local expansion history as a function of the global expansion 
factor. In all panels the solid black line represents the Hat background model, whereas the red 
dashed and blue dash-dotted lines represent an over- or underdense region. Top left: Ratio 
of the local and global Hubble rate. Top right: Local matter density parameter. Bottom 
left: Local cosmological constant density parameter. Bottom right: Local curvature density 
parameter. 

where we wrote the Hubble constant in terms of as Hl,o = lOO/i^kHis -1 Mpc -1 . 
In fig. [4] we show the time dependence of the effective local expansion history. At early times 
the curvature is negligible and the effective local Universe approaches the flat background 
Universe. At late times, the cosmological constant dominates and thus the contribution of 
matter and curvature to the energy budget becomes irrelevant. 

C Fermi Coordinates for Plane- Wave Perturbed FRW 
Universe 

Let us assume we start in Newtonian gauge with $(f G ,t 5 ) being a plane wave with wave- 
number that we can take without loss of generality in the x-direction 

®(x G ,t G ) = $ e lkXG . (109) 
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We can find the Fermi coordinates around the origin by working as in the main text and 
assume a change of coordinates valid at cubic order in the spatial distance of the form: 



t G = t L - \H(t L y L - J L $(o, t')dt' + gi (t L y L + 



(110) 
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+ fa,tr(t L ) + /io;i(^)4 + Pa^A^L^L + Z^; (^4^4 " 

This represents the most general change of coordinates around x G = at cubic order in the 
distance from the origin, and it is a straightforward generalization of (12). The subscript a 
represents the fact that those functions are zero in the limit of isotropic perturbations. By 
imposing the metric in the new local coordinates to be of the Fermi form, we can proceed and 
identify the unknown functions. We skip the details associated to straightforward algebra, 
and just quote the final result. In order to limit the size of the expressions, we simply quote 
the simple expressions that are obtained after we restrict to the case of constant £. 
Under the following change of coordinates 
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we obtain the following metric components: 
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We see that in the anisotropic case, the metric has non- vanishing Oi components at order r\. 
This form of the metric is important if we are interested in evaluating the bias for a non-scalar 
quantity, for which case the problem can not be reduced to the spherically symmetric case. 
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D Growth in Presence of a Long Mode 



In the main text, we derived a change of coordinates that is valid in a small region around 
a given time-like geodesic and that allowed us to describe the effect of a long wavelength 
fluctuation effectively as a local closed FRW Universe. This change of coordinates is valid at 
linear order in the long mode and at any order in the short wavelength perturbations. In the 
main text we focus on collapsed objects, as our main interest is extracting information about 
halo bias. Therefore we follow the short scale power well into the non-linear regime. On the 
other hand, the mapping can also be used to analytically examine the coupling of linear short 
wavelength modes to long wavelength modes while the short modes are still in the quasi-linear 
regime. In this regime, we are now going to explicitly compare results derived in our formalism 
to the ones obtained in standard perturbation theory. Since standard perturbation theory 
is performed in the Newtonian limit (k/aH ^> 1), we will adopt the simplifying assumption 
that the long mode is sufficiently far inside the horizon that the Newtonian approximation 
holds also for the long mode. In this limit we can for example neglect $ <C V 2 $/(a 2 if 2 ) 
We will also restrict ourselves to the Einstein de-Sitter (EDS) Universe for simplicity. 

As shown in the main text, the effect of a long wavelength mode on the local dynamics can 
be ascribed to a non- vanishing spatial curvature K in a fictious closed global FRW Universe. 
The curvature parameter and the curvature K are related by 

(H3) 

Let us begin to investigate the growth of short scale fluctuations in a closed FRW Universe. 
The linear growth equation for the short wavelength matter density perturbations reads as 

5 S + 2H5 S - AnGpds = . (114) 

This equation is solved by the linearly growing modes 5 s (t) = D(t)8 Bj o' 

D(t) = ln m H%H{t) / —^-3 , (115) 
2 Jo [aH(a)\ 

which in EDS simplifies to D{t) = a(t), where we use the subscript o to indicate present 
time and where we have normalized ao = 1- We will now look at the relation between 
the linear growth in a globally flat Universe and the effective local curved Universe. Using 
H(t) = —3H 2 (t)/2 in EDS, we obtain for the local effective curvature in terms of the long 
wavelength density perturbation 



K 



l^a^o . (116) 



,r G rG 



Here we have inserted a subscript i to 5 to make it more explicit that it represents a long 
wavelength fluctuation. The growth now depends on the effective curvature in two ways. First, 

12 We stress that we perform this approximation just in this appendix to make contact with former literature, 
but we do not do this same approximation in the main text, where the derivation is performed in full GR. 
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the growth in overdense regions is enhanced by a factor of 20^/21. Furthermore, from (97) we 
obtain for the relation between local and global expansion a^t/,) = (1 — 5i/3)ciG(tL), and, at 
a given fixed proper time, thus we have to evaluate the local growth at an earlier (later) scale- 
factor for overdense (underdense) regions. This partially cancels the first dependence. Adding 
both contributions, the derivative of the growth rate with respect to the long wavelength 
density reads as 



dD 



d5n 



dD 



dtt 



K 



K 



OS 



1,0 



+ 



dD 



da ] 



da j 



ao 



20 

21 ( 



1 



>'G 



13 
21 ( 



l G i 
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Thus, we finally have with 5i(t) = aaSifi 



D(6i ^0) = D C 



dD 



$l,o = a G 



21 
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It turns out that the coupling strength of 13/21 is a particlular property of the Einstein-de- 
Sitter Universe. In a more general ACDM Universe the coupling is less strong. Thus we will 
write the enhanced growth generally as D = D (l + (35i) in the following. From the rescaling 
between the global and local Hubble rate we can derive the rescaling of the local mean density 



Hj + K/a 1 



PL 
PG 



leading to 



(119) 
(120) 

Pl(x) when evaluated at the same 
physical point. Given that in this approximation we are neglecting the difference between ti 
and tG) we have 



Pl = Pg(1 + Si) 
Since p is a scalar, local and global density agree pg(x) 



5l(x) 



p(x) 



- 1 



5 G (x) 



p(x) 



G 



PL PG 

Manipulating the last expression, we obtain 

S G (x) = M 1 + + Si) +6t = 8 Li0 [l + (1 + 0)5 t ] + Si 



[l-6i)(l + 6 L )-l. 



:i2D 



>L,0 



1 34 r 

1 H Si 

21 ' 



+ Si , (122) 



where in the last step we have assumed EDS Universe. Here, we accounted both for the 
excess growth in the local frame and for the rescaling of the local mean density, with respect 
to which the local overdensity is defined. The three point function between long and short 
modes thus reads 

34 

(S G , s (x)5 G , s (x)5i(x)) = 2 x —cr 2 s Pi(k) . (123) 
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D.l Correlators between Long and Short Modes 

The coupling between long and short modes in the Newtonian regime can also be examined 
using perturbation theory (for a review see [39]). Standard perturbation theory solves the 
Newtonian fluid equations using a perturbative expansion in matter density and velocity 
divergence. In an Einstein de Sitter Universe, the second order contribution to the matter 
density field can be calculated as 

5 {2) (k) = J 0- 3 F 2 (q,k-q) 5^(q) 6^(k-q) (124) 

where 5^' is the linearly evolved primordial density field and the second order mode coupling 
kernel is defined as 



We can now apply (124) to the case where we have a Universe with short modes S s (k s ) and a 
spherical symmetric monochromatic long mode Si(ki). In this case we get for the matter field 
up to second order 

= J ^F 2 (k s , k^fyS^ik) = (l + g^)) , (126) 

where we neglected the coupling of the short and long modes with themselves. 
The skewness of the density field at second order is 

= 6 / (gs / 0ji P (l)PW)m, <?) = 3 x f^' . (127) 

The prefactor 3 arises from the fact that all three density fields in (5 3 ) can be expanded to 
second order. For the correlator between short and long modes we obtain 

(<Ux)5 s (f)M*)) = 4 1 ^ j 0^Ps(q)W)FMq) = 2 x ^lP{h) , (128) 

where we assumed Pi(q) = (27r) 3 <5( D )((f— k)P(k), k\ is the long wavelength and the prefactor 
2 arises from the fact that now only two of the three fields can be expanded to second order. 
This is in perfect agreement with our result in (123). 



E Spherical Collapse Dynamics 

The collapse of a dark matter halo can be calculated considering a spherical overdensity 
within an otherwise homogeneous background Universe. In the standard calculation the 
background is assumed to be a flat matter-only Universe (aka Einstein-de-Sitter Universe). 
After reviewing the standard spherical collapse dynamics we extend the calculation to the case 
where the background Universe is curved. As we argued in the main text, this corresponds to 
the collapse in the presence of a long-wavelength mode. This procedure will offer us a way to 
match our General Relativistic definition of the bias with the standard Newtonian definition. 
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E.l Collapse in Flat FRW 

According to Birkhoff 's theorem, a spherically symmetric overdense region evolves as a closed 
FRW Universe, whose Friedmann equation reads as 

2 _ ( a c V _ 8nGpc K c 
H °-{n) --Sr-^' (129) 

where the subscript C is used to refer to the collapsing region. This collapsing region typically 
has the size of a dark matter halo and should not be mistaken for the local patch described 
in the main part of this paper, which can contain many of these collapsing regions. The time 
evolution of the scale factor of the closed patch can be parametrized by the cycloid solution 

a c = A c (l -cos#) , t = B c (9-sm6) , with 9 E [0, 2n] , (130) 

where we defined 

3#c 3K C /2 

The spherical overdense region described by this parametrization expands until 9 = ir, then it 
turns around to collapse at 9 — 2ir, corresponding to the collapse time t co u = 2irBc- Formally 
the expansion at the collapse time is zero, but physically one expects the region to form a 
virialized object at some time between turnaround and collapse. At early times ^ < 1 the 
parametric solution can be expanded as 



a c =A c —[l + + ... , (132) 

2 V 12 360 20160 J y ' 

9 3 ( 9 2 9 4 9 e \ , , 

t =B C - 1 - — + — - 7 —- + ...) . (133) 



Solving the above equations consistently up to order 0{9 A ) one obtains 

6 2/3 / t X 2/31 



6 2/3 f t ^ 2/3 

a c = A c — 



B, 



c 



1 - 



20 \B C 



(134) 



The linear overdensity of the closed Universe collapsing at t co u is then given by the fractional 
deviation between the local and the background volume (described here by the respective 
expansion factors) 



a 



3 3 /37T \ 2/3 f_t_\ 2/3 3 (3n\ 2/3 1 + z coll 



a 6 c 5 V 2 



* t „„) = ^-l = ,|-j =l^f j^, (135, 



where we have used that the matter-only background Universe evolves according to as oc 
Finally, one obtains the critical density for collapse at z co u, linearly extrapolated to the present 
time z(t ) = 

3 /37r\ 2/3 

Screen) = 5 ( y J (! + z «xi) ~ L686 (! + Zcou) ■ (136) 
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E.2 Closed Background 

We will now extend the above calculation to the case, where the collapsing region resides 
in a curved background Universe following [H] and [45] . We will consider the case of an 
overdense, closed background Universe and note that the open background can be treated 
analogously. Furthermore, we will restrict ourselves to a background Universe without a dark 
energy component, such that only matter and curvature contribute to the energy budget. 
This closed background Universe is parametrized as 



a B 



A 



B 



COS if) 



t = B B (tj — sin 77) 



with i] G [0, 2tt] and the parameters 

AkGpb 



A 



B 



3K 



and B 



B 



AirGp 



B 



B 



3K 



3/2 



(137) 
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This curved background can now be identified with the effective curved patch describing a 
long wavelength fluctuation. We study the evolution of a collapsing spherical overdensity, 



expanding the parametric solutions for both the background (137) and the collapsing re- 



gion (130) at early times. This means that we restrict ourselves to treat the curvature of the 



background at linear order. The linear density contrast then scales as 



- 1 



3 /3tt\ 
5 IT J 



2/3 



1 

*coll 



2/3 



2/3' 



t 2/3 



(139) 



where to = 2nB. For K — > we have tn — > 00 and thus we recover the EDS result shown 
above. The overdensity for an object that collapses at t co u, linearly extrapolated to present 
time thus reads as 



[Zcoll 



1 + Zcoll 



3tt 

Y 



2/3 



tn) 



2/3 



tcoll \ 



2/3" 



We can now write down the collapse time for the background Universe. 

TlQ m B 7T (1 — VL K n) 



t n = 2nB B 



(140) 



(141) 



Using (toAcoii) 2/3 
rewritten as 



1 + Zcoii and to ~ 2/ (3H ) the overdensity of the collapsing region can be 



S c {z co ii) _ 3 / 3_7r\ 
1 + z coll ~ 5 V 2 J 



2/3 



1 



3tt 



2/3 



n 



K,B 



1 



;i - n K>B y/z 1 + z co u 



(142) 



We are now going to consider the case where the curvature of the background Universe can 
be described by a long wavelength fluctuation with present day amplitude S^b- In this case 
we obtain for the density parameters of the background Universe 



n 



1 + 5, 



LB 



m,B 



(1 - <W3) 2 



LB 



K.B 



76, 



LB 



(143) 
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For the overdensity of a perturbation that collapses at z coa linearly extrapolated to the present 
day we obtain 

5 c (z coll ) « 1.686 (1 + z coll ) - 5 lfi = 5 c (n K>B = 0) - 8 hB ■ (144) 

This result is of course very intuitive and it allows us to explicitly verify that our General 
Relativistic definition of the bias agrees with the standard Newtonian one. 



F Perturbed Geodesic Parameters 

In this appendix we will provide the essence of cosmological perturbation theory and expalain 
the gauge choices used in this paper. Then we will quickly review the most important formulae 
required for the mapping to observables before we conclude by specialising our result for the 
observed overdensity to the case where the matter distribution itself is the tracer. 

F.l Gauge Transformations 

The most general perturbed metric for a flat Universe reads as [UJ 

ds 2 = -(1 + 2A)dt 2 - 2aB i dx i dt + a 2 [(1 + 2L% i + E tj ] dx i dx j . (145) 

We will restrict ourselves to scalar modes Bi = BQf^ and = EQ^\ where is the 
scalar eigenmode of the Laplacian. Here we perform the scalar-vector-tensor decomposition 
in fc-space, where = exp ik • x . We consider a Universe filled with dark matter plus 
dark energy and neglect anisotropic stress and pressure perturbations. 

A gauge transformation corresponds to a change in spatial position x % and comoving time 
adr = dt 

x* = x { + V f = T + T (146) 

under such a transformation the metric perturbations transform as 

k 

A = A-aT- aHT D = D - -L - aHT (147) 

B = B + aL + kT = + kT 

and the components of the energy momentum tensor transform as 

5 = 5 + 3aHT v = v + at , (148) 

where Vi = vQ\ . Gauge invariance refers to the fact that certain combinations of metric 



and energy momentum perturbations are invariant under a change of coordinates (146), i.e., 
the numerical value of a quantity does not change. Gauge invariance is necessary, but not 
sufficient, for observability. 
We will consider two gauges 
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1. Newtonian Gauge 
Newtonian gauge is denned by B = E = setting A = $ and D = . 

ds 2 = -(l + 2$)rft 2 + a 2 (l-2^)dx 2 (149) 

Neglecting anisotropic stress we have $ = \P. For the Einstein equations we have 

-k 2 <$> - 3a 2 H 2 [ $ + ^ J =47rGa 2 /^ (150) 

aH\$+— ) =4nGa 2 (p + p) — = -a 2 H— (151) 

2. Comoving Gauge 

In comoving gauge we set E = and 5Tf = corresponding to Vi = B^. Setting A = £ 
and D = ( leads to the metric 

cfe 2 = -(1 + 2£)cft 2 + avidx'dt + a 2 (t)(l + 2C)df 2 (152) 

then the Einstein equations read as 

/ ? .(com)\ 

A; 2 f ( + aH —^~ ) =4vrGa 2 p m 5 {com) (153) 

- C =0 , (154) 

where t>( com ) = v^ N \ One can show that on scales larger than the sound horizon £ is 
constant. Using £ = we see that the lapse function £ vanishes in pressureless media 
and thus the comoving gauge is also synchronous, i.e., proper time agrees with the 
coordinate time. Using 

V (N) _ B (N) V (N) 
Tn-^cotu = 7 = 7 (155) 

the overdensities in the Newtonian and comoving gauge are related by 

5 {com) = § (N) _ a p TN _ ¥com = S W + 3aR ( 156 ) 

k 

For the spatial metric perturbations we have 

V W H 2 ( $\ 

C = - aHT N ^ corn = -$ - a# — = -$ + U + - j (157) 

In synchronous gauge A = 5 = there are no sources in the equation of motion for 
the velocity of stress free matter, i.e., if it was at rest initially it will remain so for all 
times |48] . In this case the density perturbation in comoving and synchronous gauge 
agree <5( com ) = §( s v n )_ The CMBFAST Boltzmann code (23] is providing the synchronous 
gauge transfer function and can thus be used to infer the transfer function for the 
comoving gauge density perturbation. 

Combining the Einstein equations in comoving and Newtonian gauge we have 

- k 2 <& = 4vrGa 2 (5 (com) , (158) 

which is valid on all scales. 
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F.2 Volume Distortion Sz. Observed Redshifts 



Here we explain the symbols used in eq. (50) for the reader's convenience (see pQ for a detailed 
explanation and derivation): 



J = SV/V = - $ + v% - (1 + z)-^6zg^ z - 2^5z G ^ z - 5z G ^ z 



1 + z dH r 5r 
-2k + ——--5zg^ z + 2— 
H dz r 



(159) 
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8 JO 



dlnH 
dln(l + z) 



dr^-2 



1 + z 
Hr« 



v l e; - $ - 2 



dr a<3? 



T — T * 1 

dr — V$ + — 

rr s H 



1 dvie 1 
a dr 



In the evaluation of the above expression we used that the total derivative is given by d/dz 
H~ x d/dr = —(d — e l di) = — (d — d r ) and that the velocity follows the evolution equation 



Vi + Hvi = -- 
a 

The perturbation to the redshift of the source 5zg^ z , is given by the relationship 

5z G ^ z = aH5r + [vie 1 — <&]* — 2 / a<j> dr . 



(160) 



(161) 



Here the four velocity of the source is given by u a = a _1 ((l + $),?/) and e l is the photon 
propagation direction as seen from the observer, r is the comoving line-of-sight distance, r s 
is the comoving line-of-sight distance of the source. 5t q is the perturbation to the conformal 
time at the time of observation, and it is just a monopole term that cancels when measuring 
fluctuations. 5r is the radial displacement, given by 



5r = 8tc\ 



dr $ 



(162) 



The deflection of the photons on their way form the source galaxy to the observer can be 
quantified as 



59 

8(f) 



dr 
dr 



r s — r 

r r s 

r s — r 
rr, sin^ 



(163) 



The latter can be combined to calculate the distortion of the solid angle as quantified by the 
convergence k 



k 



dr 



r s — r 



V 2 $ 



(164) 



where V is the differential operator on the two dimensional unit sphere. In case the survey 
is not volume limited but rather flux limited, we need to replace J with 



J J-hp5Vi 



(165) 
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where 8T>l is the perturbation to the luminosity distance, which is given by 



SV L (z) 

Dr(z) 



1 + Vi e l - - ^-^<f*G_+* +(n o +-) St,, 



+2 / dr 
Jo 



$ r • (r, — r)r /„, d(a$) 

a$ + 1 * ; V 2 $ - a^— -L + 2a*,; e l 

r„ r„ 2r \ at 



(166) 



where D L (z) = (1 + z)r(z) is the unperturbed luminosity distance and p is the slope of the 
luminosity function. 

Using dz = —H(l + z)dt we can write in a general and in Newtonian gauge 



5z G . 



^g(^l) - Z G (tg) 

1 + zg(*g) 

*t<3 



-H(t G )(t L -t G ) = -H 1* ' dtA 

Jo 
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dt $ = — aif 
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(167) 



( ^ = iz^ =[( „_ J%)e ._ x]: 

/■r s 

= [vie 1 — <&]* — 2 / rfra$. 
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Under a gauge transformation ( 146 ) we have 



5z G ^ L = 5z G ^ L + aHT 



5z 



8zg-+z + aHT . 



(169) 



Thus Sz g ^l — &z G ^ z does not change under the gauge transformation and is gauge invariant. 



For the evaluation of the full expression in eq. (61) we first transform all the quantities to 



fc-space. First, we consider the line of sight projection of the velocity 



v 4 e 



d 3 k 
(2^)3 



(-i^ N \k)) Q^(jfc) 



(0)/ 



(170) 



where \x = x ■ k/(xk) is the cosine between the fc-mode and the line of sight. For the redshift 
space distortion term we have then 



d r Vie % = n 3 dje l Vi 



Thus the volume distortion term reads as 



d 3 k 
(2tt) 3 



{^ N \k)) Q(°>(£) 



(171) 



J(k) = -$(£) + A(z)(-i^ N \k) - $(£)) + (f(z) - - V 2 ^v (N) (k) . (172) 

The full observed density perturbation is the sum of the latter and the perturbation in the 
proper number density of tracers 

S p (k) = b5 {com) (k) + b c ({k) + B(z) (-^v (N \k) - ifiv {N \k) - $(jfe) J , (173) 
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P(z) aH 



(174) 



where we related the density and velocity perturbations to the Newtonian gauge metric per- 
turbation 



5 {com) (k) 
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aH J AnGp 
and introduced the auxiliary functions 
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d log n p 



iJ 2 W 01og(l + *) ' 

The power spectrum is then given by (2n) 3 5^ D \k + k')P obs (k) = (^S b s (k)5l bs (k')^ . 
F.3 Matter as Tracer 

For matter we have n p (t L ]Vt K ) = p(t L )(l + D(t L )5^ om ^) . This can be seen in two ways: 
firstly, in synchronous slicings the proper time and the coordinate time agree, thus the matter 
overdensity in the local frame must agree with the matter overdensity in comoving gauge. 



Also, we saw above in eq. (102), that the local matter density is related to the global one by 



Pl = Pg(1 + D{ti)5 com fl)- Using that VLk,o oc 5 com ,o and that <91ogp/<91og(l + z) = 3 we have 

5 obs (z, 9, 0) = <f (com) (s G ) + 3fc M - 35z G ^ z + J , (177) 
where the Sz's are in a general, yet unspecified gauge. 

For the transformation from a comoving to general gauge we have for the densities 



(178) 



6 ( 9 en) = S (com) + 3 a JJT com ^ gen . 

The integral entering into Szq^l transforms as 

38z G 1 ™ ) L = -3H J dtA {9en) = -3H J dtA {com) +3aHT c 

(179) 

where we used that A^ 00 ™) = and solved eq. (178) for T com ^ gen . Evaluating 5 b s in the 
general gauge we obtain 



— 3fl HT — X(9en) _s:(com) 



S obs (z, 0, <p) = 5 {com) - 3H J dtA^ - 35 z, 



{gen) 
G^tz 



J = «f&«0 - 38z^l + J, 



;iso) 
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where S^ 9en ^ and 5zq™[ are the matter overdensity and the redshift lapse in a general gauge. 
This expression agrees with 5 np = 5^ - 35 z^ in pQ. In the follow up paper, [13] used a 
local bias in the matter density at the observed redshift 5^ 9en ^ — 35zq™ z . This means that in 
contrast to our approach the bias factor b is also multplying the redshift lapse terms 
between global, local and observed redshift, i.e. the evolution of the sample is fixed to be 
<91ogn p /<91og(l + z) = 3b, while for a typical quasar sample this number can vary in a much 
wider range depending on the redshift distribution of the sample. 



References 

[1] J. Yoo, A. L. Fitzpatrick and M. Zaldarriaga, "A New Perspective on Galaxy 

Clustering as a Cosmological Probe: General Relativistic Effects," Phys. Rev. D 80 
(2009) 083514 [arXiv:0907.0707 [astro-ph.CO]]. 

[2] N. Dalai, O. Dore, D. Huterer and A. Shirokov, "The imprints of primordial 

non-gaussianities on large-scale structure: scale dependent bias and abundance of 
virialized objects," Phys. Rev. D 77 (2008) 123514 [arXiv:0710.4560 [astro-ph]]. 

[3] A. Slosar, C. Hirata, U. Seljak, S. Ho and N. Padmanabhan, "Constraints on local 
primordial non-Gaussianity from large scale structure," JCAP 0808 (2008) 031 
[arXiv:0805.3580 [astro-ph]]. 

[4] L. Senatore and M. Zaldarriaga, "The Effective Field Theory of Multifield Inflation," 
arXiv:1009.2093 [hep-th]. 

[5] T. Baldauf, U. Seljak and L. Senatore, "Primordial non-Gaussianity in the Bispectrum 
of the Halo Density Field," JCAP 1104 (2011) 006 [arXiv:1011.1513 [astro-ph.CO]]. 

[6] E. Fermi, "Sopra i fenomeni che avvengono in vicinanza di una linea oraria," Atti 
Accad. Naz. Lincei CI. Sci. Fis. Mat. Nat. 31,184-187, 306-309, 1922. F. K. Manasse 
and C. W. Misner, "Fermi Normal Coordinates and Some Basic Concepts in 
Differential Geometry," J. Math. Phys. 4, 735 (1963). 

[7] G. Tormen and E. Bertschinger, "Adding Long Wavelength Modes to an iV-Body 
S imulat ion , " arXiv : astro-ph /9512131. 

[8] S. Cole, "Adding Long- Wavelength Power to N-body Simulations," 
arXiv:astro-ph/9604046. 

[9] M. D. Schneider, S. Cole, C. S. Frenk and I. Szapudi, "Fast generation of ensembles of 
cosmological N-body simulations via mode-resampling," arXiv: 1103.2767 [astro-ph.CO]. 

[10] M. Sasaki, "The Magnitude - Redshift relation in a perturbed Friedmann Universe," 
Mon. Not. Roy. Astron. Soc. 228, 653-669 (1987). 



45 



[11] S. Dodelson, F. Schmidt and A. Vallinotto, "Universal Weak Lensing Distortion of 

Cosmological Correlation Functions," Phys. Rev. D 78, 043508 (2008) [arXiv:0806.0331 
[astro-ph]]. 

[12] D. Wands and A. Slosar, "Scale-dependent bias from primordial non-Gaussianity in 
general relativity," Phys. Rev. D 79 (2009) 123507 [arXiv:0902.1084 [astro-ph.CO]]. 

[13] J. Yoo, "General Relativistic Description of the Observed Galaxy Power Spectrum: Do 
We Understand What We Measure?," Phys. Rev. D82, 083508 (2010). 
[arXiv: 1009.3021 [astro-ph.CO]]. 

[14] N. Bartolo, S. Matarrese, A. Riotto, "Relativistic effects and primordial 

non-Gaussianity in the galaxy bias," JCAP 1104, 011 (2011). [arXiv: 1011.4374 
[astro-ph.CO]]. 

[15] E. Komatsu et al, "Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) 
Observations: Cosmological Interpretation," arXiv: 1001. 4538 [astro-ph.CO]. 

[16] L. Senatore, K. M. Smith and M. Zaldarriaga, "Non-Gaussianities in Single Field 

Inflation and their Optimal Limits from the WMAP 5-year Data," JCAP 1001 (2010) 
028 [arXiv:0905.3746 [astro-ph.CO]]. 

[17] P. Creminelli, M. A. Luty, A. Nicolis and L. Senatore, "Starting the Universe: Stable 
violation of the null energy condition and non-standard cosmologies," JHEP 0612 
(2006) 080 [arXiv:hep-th/0606090]. P. Creminelli, G. D'Amico, J. Norena and 
F. Vernizzi, "The Effective Theory of Quintessence: the w < — 1 Side Unveiled," JCAP 
0902 (2009) 018 [arXiv:081 1.0827 [astro-ph]]. P. Creminelli, G. D'Amico, J. Norena, 
L. Senatore and F. Vernizzi, "Spherical collapse in quintessence models with zero speed 
of sound," JCAP 1003 (2010) 027 [arXiv:091 1.2701 [astro-ph.CO]]. 

[18] N. E. Chisari and M. Zaldarriaga, "On the connection between Newtonian simulations 
and General Relativity," arXiv: 1101.3555 [astro-ph.CO]. 

[19] A. L. Fitzpatrick, L. Senatore and M. Zaldarriaga, "Contributions to the Dark Matter 
3-Pt Function from the Radiation Era," JCAP 1005 (2010) 004 [arXiv:0902.2814 
[astro-ph.CO]]. 

[20] C. Cheung, P. Creminelli, A. L. Fitzpatrick, J. Kaplan and L. Senatore, "The Effective 
Field Theory of Inflation," JHEP 0803 (2008) 014 [arXiv:0709.0293 [hep-th]]. 

[21] G. Tormen and E. Bertschinger, "Adding Long Wavelength Modes to an A r -Body 
Simulation," arXiv:astro-ph/9512131. 

[22] S. Cole, "Adding Long- Wavelength Power to N-body Simulations," 
arXiv:astro-ph/9604046. 



46 



[23] U. Seljak and M. Zaldarriaga, "A Line of Sight Approach to Cosmic Microwave 
Background Anisotropics," Astrophys. J. 469 (1996) 437 [arXiv:astro-ph/9603033]. 

[24] A. Lewis, A. Challinor and A. Lasenby, "Efficient Computation of CMB anisotropies in 
closed FRW models," Astrophys. J. 538 (2000) 473 [arXiv:astro-ph/9911177]. 

[25] A. Gangui, F. Lucchin, S. Matarrese and S. Mollerach, "The Three Point Correlation 
Function Of The Cosmic Microwave Background In Inflationary Models," Astrophys. J. 
430 (1994) 447 [arXiv:astro-ph/9312033]. 

[26] D. H. Lyth, C. Ungarelli and D. Wands, "The primordial density perturbation in the 
curvaton scenario," Phys. Rev. D 67, 023503 (2003) [arXiv:astro-ph/0208055]. 

[27] M. Zaldarriaga, "Non-Gaussianities in models with a varying inflaton decay rate," 
Phys. Rev. D 69, 043508 (2004) [arXiv:astro-ph/0306006]. 

[28] P. Creminelli and L. Senatore, "A smooth bouncing cosmology with scale invariant 
spectrum," JCAP 0711 (2007) 010 [arXiv:hep-th/0702165]. 

[29] X. Chen and Y. Wang, "Quasi-Single Field Inflation and Non-Gaussianities," JCAP 
1004 (2010) 027 [arXiv:091 1.3380 [hep-th]]. 

[30] M. Bruni, R. Crittenden, K. Koyama, R. Maartens, C. Pitrou and D. Wands, 
"Disentangling non-Gaussianity, bias and GR effects in the galaxy distribution," 
arXiv: 1106.3999 [astro-ph.CO]. 

[31] E. Silverstein and A. Westphal, "Monodromy in the CMB: Gravity Waves and String 
Inflation," Phys. Rev. D 78 (2008) 106003 [arXiv:0803.3085 [hep-th]]. 

[32] L. McAllister, E. Silverstein and A. Westphal, "Gravity Waves and Linear Inflation 
from Axion Monodromy," arXiv:0808.0706 [hep-th]. 

[33] D. Green, B. Horn, L. Senatore and E. Silverstein, "Trapped Inflation," Phys. Rev. D 
80 (2009) 063533 [arXiv:0902.1006 [hep-th]]. 

[34] N. Barnaby, Z. Huang, L. Kofman and D. Pogosyan, "Cosmological Fluctuations from 
Infra-Red Cascading During Inflation," Phys. Rev. D 80 (2009) 043501 
[arXiv:0902.0615 [hep-th]]. 

[35] R. Flauger, L. McAllister, E. Pajer, A. Westphal and G. Xu, "Oscillations in the CMB 
from Axion Monodromy Inflation," arXiv:0907.2916 [hep-th]. R. Flauger and E. Pajer, 
"Resonant Non-Gaussianity," arXiv: 1002.0833 [hep-th]. 

[36] X. Chen, R. Easther and E. A. Lim, "Large Non-Gaussianities in Single Field 
Inflation," JCAP 0706 (2007) 023 [arXiv:astro-ph/06 11645]. 



47 



[37] C. M. Hirata and U. Seljak, "Analyzing weak lensing of the cosmic microwave 
background using the likelihood function," Phys. Rev. D 67 (2003) 043001 
[arXiv:astro-ph/0209489] . 

[38] C. M. Hirata and U. Seljak, "Reconstruction of lensing from the cosmic microwave 
background polarization," Phys. Rev. D 68 (2003) 083002 [arXiv:astro-ph/0306354]. 

[39] F. Bernardeau, S. Colombi, E. Gaztanaga and R. Scoccimarro, "Large-scale structure 
of the Universe and cosmological perturbation theory," Phys. Rept. 367 (2002) 1 
[arXiv:astro-ph/01 12551] . 

[40] A. Challinor and A. Lewis, "The linear power spectrum of observed source number 
counts," arXiv: 1105.5292 [astro-ph.CO]. 

[41] C. Bonvin and R. Durrer, "What galaxy surveys really measure," arXiv:1105.5280 
[astro-ph.CO]. 

[42] D. Jeong, F. Schmidt, C. M. Hirata, "Large-scale clustering of galaxies in general 
relativity," [arXiv: 1107.5427 [astro-ph.CO]]. 

[43] D. Seery and J. E. Lidsey, "Primordial non-Gaussianities in single field inflation," 
JCAP 0506 (2005) 003 [arXiv:astro-ph/0503692]. X. Chen, M. x. Huang, S. Kachru 
and G. Shiu, "Observational signatures and non-Gaussianities of general single field 
inflation," JCAP 0701 (2007) 002 [arXiv:hep-th/0605045]. L. Boubekeur, P. Creminelli, 
J. Norena and F. Vernizzi, "Action approach to cosmological perturbations: the 2nd 
order metric in matter dominance," JCAP 0808 (2008) 028 [arXiv:0806.1016 
[astro-ph]]. 

[44] Peebles, P. J. E. 1980. The large-scale structure of the Universe. Princeton University 
Press, 1980. 435 p. . 

[45] C. G. Lacey and S. Cole, "Merger rates in hierarchical models of galaxy formation," 
Mon. Not. Roy. Astron. Soc. 262, 627 (1993). 

[46] F. I. Cooperstock, V. Faraoni and D. N. Vollick, "The influence of the cosmological 
expansion on local systems," Astrophys. J. 503, 61 (1998) [arXiv:astro-ph/9803097]. 

[47] W. Hu, "Lecture Notes on CMB Theory: From Nucleosynthesis to Recombination," 
arXiv:0802.3688 [astro-ph]. 

[48] C. P. Ma and E. Bertschinger, "Cosmological perturbation theory in the synchronous 
and conformal Newtonian gauges," Astrophys. J. 455, 7 (1995) 
[arXiv:astro-ph/9506072] . 



48 



